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PREFÁCIO 


O presente livro é o primeiro volume de uma série intitulada “Práticas em Saneamento, 
Clima e Meio Ambiente” publicado pela Editora Poisson, apresentando métodos práticos 
para aulas em campo e laboratório considerando a importância do ensino aplicado em 
engenharia e áreas afins. 


O objetivo central deste livro foi apresentar a aplicação de diversos aspectos 
metodológicos para obtenção de dados primários em campo e laboratório, os quais são 
úteis para subsidiar acadêmicos, pesquisadores, engenheiros e técnicos da Região 
Amazônica. 


Assim, em razão da escassez de informações práticas que se tem nessa região e sabendo 
o quão importante é a divulgação científica, evidencia-se a necessidade de promover, na 
forma de coletânea técnica, a elaboração de dez capítulos, os quais são descritos a seguir: 


Capítulo 1. Aula prática: medição de vazão; 

Capítulo 2. Aula prática: medição de vazão usando drones; 
Capítulo 3. Aula prática: infiltração da água no solo; 
Capítulo 4. Aula prática: hidrologia estatística; 


Capítulo 5. Aula prática: detecção de mudanças a partir de série temporal de imagens 
orbitais utilizando os métodos diferença de NDVI e RCEN modificada; 


Capítulo 6. Aula prática: nivelamento geométrico, determinação de desníveis e altitudes; 
Capítulo 7. Aula prática: sistema de recalque levantamento da curva de uma bomba 
centrífuga; 

Capítulo 8. Aula prática: carga fixa-poluição; 

Capítulo 9. Aula prática: ensaio de adensamento por gravidade em batelada de lodo de 
ETA; 


Capítulo 10. Aula prática: determinação de oxigênio dissolvido em águas, método 
químico. 


A disseminação dessa coletânea técnico acadêmica, realizada arduamente com as 
experiências práticas e acadêmicas dos diversos professores que contribuíram em cada 
capítulo, desenvolvida de maneira detalhada e didática, para que possa ser replicada em 
outras instituições de ensino, bem como ser usada como ferramenta para obtenção de 
dados de práticos, propiciando assim à transferência dos conhecimentos e resultados 
não somente na Região Norte, mas em todo o Brasil. 
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Capítulo 1 
Aula prática: Medição 


de vazão 


Giovanni Chaves Penner 
Katiucia Nascimento Adam 
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1.1. OBJETIVOS INSTRUCIONAIS 
Ao final desta aula prática de medição de vazão em pequenos cursos d'água você deverá ser capaz de: 
e Enunciar três possíveis métodos aplicáveis a medição de vazão em pequenos cursos d'água; 
e Conceituar batimetria, molinete e traçador; 
e Operacionalizar a medição de vazão um pequeno curso d'água; 
e Identificar as incertezas do método de medição de vazão empregado; 


e Comparar os resultados obtidos por diferentes métodos de medição de vazão identificando 
resultados anômalos. 


1.2. MEDIÇÃO DE VAZÃO EM HIDROLOGIA 


A visão de engenharia dada a nossa disciplina de Hidrologia e Climatologia foca no o ciclo hidrológico do 
ponto de vista quantitativo. Uma informação deveras solicitada em várias aplicações práticas em 
Hidrologia é “qual a vazão”. Outra dúvida frequente associada ao tema: existe variação significativa da 
vazão ao longo do ano que comprometa um dado uso? Para responder, muitas vezes, faz-se necessário 
medir a vazão do curso d'água de interesse, criando-se séries históricas que são extremamente úteis para 


diversos estudos e projetos de Engenharia. 


Desta forma mede-se vazão para se ter valores de referência para a tomada de decisão. Em muitos casos 
uma única medida não suficiente, daí a necessidade de séries históricas diárias, mensais e anuais. 
Eventualmente, em condições com as da cidade de Belém-PA, são necessárias medições horárias devido 
aos efeitos da maré. 


1.3. MÉTODOS DE MEDIÇÃO DE VAZÃO EM PEQUENOS CURSOS D'ÁGUA 


Como foi apresentado em sala de aula há diversos métodos para medir vazão em cursos d'água, cada um 
com condições de aplicação e ferramental específico. Como o foco da nossa aula prática é para pequenos 
cursos d'água serão considerados os métodos descritos a seguir. 


1.3.1. FLUTUADOR 


O método consiste em determinar a velocidade de deslocamento de um objeto flutuante, medindo o tempo 
utilizado para o seu deslocamento num determinado trecho de rio de comprimento conhecido (SANTOS et 
al., 2001). O flutuador de superfície pode ser uma boia, uma fruta, uma garrafa plástica, ou qualquer outro 
dispositivo. Um flutuador de superfície move-se com a mesma velocidade da superfície da água, e se for 
muito leve pode sofrer a influência do vento. Por outro lado, como o flutuador indica apenas a velocidade 
da superfície da água, é necessário aplicar um coeficiente redutor para obter a velocidade média na seção 
(SANTOS et al., 2001). A literatura (AZEVEDO NETTO et al., 1998) recomenda que a velocidade média seja 
obtida multiplicando-se a velocidade superficial por um coeficiente entre 0,8 e 0,9 (0,8 para rios com 
fundo pedregoso ou 0,9 para rios com fundo barrento). Todavia existem casos que indicam valores 
distintos e essa pode ser uma das conclusões desta aula prática. Com a velocidade média determinada faz- 
se a batimetria da seção para que seja determinada a área de escoamento da água. Pelo produto da 
velocidade média pela área da seção de escoamento tem-se a vazão. 


1.3.2. MOLINETE 


Mais um método que determina a vazão pelo produto da velocidade do escoamento pela área da seção 
transversal. 


Molinete são pequenas hélices que giram impulsionados pela passagem da água. Os molinetes são 
instrumentos projetados para girar em velocidades diferentes de acordo com a velocidade da água. 
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Em cursos d'água naturais, a velocidade da água, em geral, é maior no centro de um rio do que junto às 
margens. Adicionalmente, a velocidade da água é mais baixa junto ao fundo do rio do que junto à 
superfície. Em função desta variação da velocidade nos diferentes pontos da seção transversal, utilizar 


apenas uma medição de velocidade pode resultar em uma estimativa errada da velocidade média. 


Para obter uma boa estimativa da velocidade média é necessário medir em várias verticais, e em vários 
pontos ao longo das verticais conforme descrito nos slides da aula e na apostila de apoio. Lembrando que 
cada vertical de medição de velocidade da água terá uma velocidade média e estará associada a uma área 
de escoamento. A soma das vazões associada a cada vertical fornecerá a vazão total. 


1.3.3. DILUIÇÃO QUÍMICA OU TRAÇADOR SALINO 


Segundo Silva et al., (2009) traçador é qualquer substância, ou partícula /entidade (química ou biológica), 
que pode ser usada para seguir, quer pontualmente ou de forma contínua, o comportamento de um 
determinado sistema ou de um componente, tal como velocidade da água ou vazão. 


Pelo exposto a medição da vazão pelo método de diluição usa um traçador. Tal técnica se aplica a 
pequenos cursos da água, porém, a velocidade e a irregularidade do leito devem ser suficientes para sua 
mistura até o ponto de determinação da concentração. Aplica-se também onde não é possível instalar 
calhas e vertedores devido ao custo ou acessibilidade, como presença de cachoeiras. Consiste em lançar no 
curso em estudo uma vazão de uma solução que não seja encontrada em grande quantidade nessas águas 
e, a jusante, medir a concentração dessa solução comparando-a com a concentração natural do curso 
(COLLISCHONN et al., 2015; PORTO et al., 2001). 


Os produtos a serem empregados como traçadores devem apresentar certas qualidades essenciais como: 
elevada solubilidade na água, presença nula ou apenas traços em águas naturais, ser atóxico e não 
corrosivo para não causar danos ao ambiente, equipamentos e nem as pessoas que o manipulam. Sua 
concentração não deve se alterar sob influência da luz, bactérias, matéria orgânica, sedimentos em 
suspensão ou outras substancias solúveis presentes no curso d'água. Não é aconselhável o uso de 
elementos químicos radioativos como traçadores, que embora apresentem maior precisão, causam danos 


ao ambiente (HINDI et al., 1998). 


Dentre os produtos mais utilizados como traçadores tem-se o Cloreto de Sódio (NaCl]), Cloreto de Potássio 
(Nak), Nitrito de Sódio (NaNO,), Sulfato de Manganês (MnSO4), Rodamina (C,oH,1CIOsN») e Bicromato de 
Sódio (NazCr,H;0). O cloreto de sódio, popularmente conhecido como sal de cozinha refinado, além de ser 
inofensivo à ecologia, também possui a propriedade de aumentar a condutância específica da água, numa 
razão proporcional a concentração do mesmo. Tal fato permite medir a condutividade elétrica ao invés de 
se medir a concentração do sal, o que facilita e agiliza o trabalho de medição da vazão (HINDI et al., 1998). 


Para a aplicação deste método existem duas condições fundamentais citadas por Hindi et al. (1998): 


1) A homogeneidade do sistema água-traçador deve ser constante para qualquer ponto da seção de 
medição: em situação de escoamento turbulento esta homogeneidade se dá em distância muito mais curta 
quando comparadas a aquelas necessárias para rios com escoamento lento e tranquilo. 


2) Conservação da matéria: não deverá ocorrer perda de massa do traçador entre os pontos de 
lançamento e de medição. Para isso é necessário determinar uma distância mínima (L), também chamada 
na literatura especializada de comprimento de mistura, entre estes dois pontos que seja suficiente para 
garantir a homogeneidade da solução e satisfazer a condição anterior. Esta distância depende da 
velocidade da água no rio, da turbulência do escoamento, da presença de obstáculos e do modo de injeção 
do sal. 


Qualquer que seja o traçador empregado existe dois processos distintos para medir a vazão por este 
método, o processo por injeção a vazão constante e por volume instantâneo (integração). A principal 
diferença entre eles é que o primeiro processo se trata da injeção do traçador a uma vazão constante, 
enquanto que no segundo, a injeção é instantânea. A vantagem do processo por integração é não precisar 
controlar a vazão de injeção da solução no rio, o que torna mais fácil a aplicação do método (KILPATRICK 
et al., 1985). Moore (2005) indica que o processo por injeção constante seja aplicado em pequenos 
córregos com vazão menor do que 100 L/s e para maiores vazões, de até 10.000 L/s, seja aplicado o 
processo de injeção instantânea. 
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1.4. PRÁTICA EM CAMPO 


Nossa aula prática acontecerá em campo fazendo medições reais, por diferentes métodos e com 
possibilidade de comparação dos métodos e resultados. 


1.4.1. SELEÇÃO DO LOCAL 


O local escolhido deve atender aos critérios apresentados em sala de aula. De forma resumida tem-se: o 
trecho escolhido para medição da vazão deve ser retilíneo (sem curvas), ter no mínimo 15 cm de 
profundidade e não ser uma área de águas lentas. Corredeiras desobstruídas são ideais. O comprimento 
deste trecho deve algo entre 10 e 20 metros para que as medições sejam realizadas com facilidade, todavia 
em casos reais deve-se trabalhar com o trecho mais retilíneo e uniforme disponível. 


Para facilitar nossa atividade foi escolhido o canal que passa na lateral do prédio da Faculdade de 
Engenharia Naval, no Setor Saúde, da Cidade Universitária Prof. José da Silveira Netto da Universidade 
Federal do Pará em Belém. 


Nas Fotos 1 e 2 é indicada a localização do canal para a realização da aula prática. 


Foto 1 - Identificação do canal de drenagem na lateral do prédio da Faculdade de Engenharia Naval. 
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1.4.2. ANTES DA MEDIÇÃO 


Separar os seguintes materiais antes de se deslocar para o local de medição: 


1) Protetor solar; 

2) Tenda ou guarda sol e cadeiras; 

3) Trenas de 3 e 50 metros; 

4) Cronômetro; 

5) Planilha de campo (Apêndice 1), prancheta, caneta e lapiseira; 

6) Sensor de Condutividade Elétrica e Molinete com sistema de aquisição de dados LabQuest 2, 


marca Vernier; 


7) Balde graduado; 


8) Frasco de Mariotte; 

9) 5 quilos de sal de cozinha (cloreto de sódio); 
10) Garrafa e/ou boia flutuadora; 

11) Macacão de borracha. 


Não esqueça de estar vestido adequadamente de calça, botas, camisa confortável, óculos de proteção solar 
e chapéu/boné. Também fiquem atentos para fazer registros fotográficos em todas as etapas. 


1.4.3. ATIVIDADES EM CAMPO 


Em grupos de 5 alunos nesta atividade de campo serão testados os três métodos de medição de vazão 
citados. Inicialmente um grupo fará a primeira medição usando o primeiro método, em seguida outro 
grupo fará a medição usando o segundo método e assim por diante. Caso haja tempo disponível todos os 
grupos farão a medição dos seus próprios dados, caso contrário os dados medidos por cada grupo 
deveram ser compartilhados com os demais grupos para a elaboração do seu respectivo relatório. Os 
métodos usados em campo seguirão a seguinte ordem: 


Primeiro, usando o flutuador superficial, testando diferentes flutuadores (garrafas e bola de isopor). 
Comece medindo com a trena a distância entre o ponto onde será lançado o flutuador e o ponto onde será 
registrado o tempo. Um aluno será o responsável pelo lançamento e outro pelo registro do tempo, o 
procedimento deve ser repetido três vezes. Calcule a velocidade superficial para cada lançamento feito. 
Faça a batimetria da seção de escoamento da água. 


Segundo, usando o molinete, conforme a largura e profundidade da lâmina d'água defina em quantas 
verticais e quantas profundidas diferentes em cada vertical precisa ser medida a velocidade de 
escoamento da água. Faça as medições das velocidades. Determine as dimensões da seção de escoamento 
da água para ser calculada a área referente a cada vertical. 


Por fim, o terceiro, usando a diluição química com traçador salino. A primeira tentativa será usando um 
pulso, isto é, tomando aproximadamente 18 L de água do próprio curso d'água acrescente um 1 kg de sal, 
misture bem até não ter sal precipitado no fundo do balde. Meça a condutividade elétrica da solução 
contida no balde. Um aluno irá derramar a solução salina num ponto a montante, e outro fará os registros 
contínuos de condutividade elétrica, identificando o pulso das leituras, a uma distância que seja suficiente 
para que ocorra a mistura completa da solução salina. Configure o Labquest 2 para fazer registros a cada 2 
segundos e defina o tempo total do experimento usando a informação do ensaio com o flutuador. 


O método da diluição química pode ser repetido seguindo parte do procedimento acima, todavia usando a 
vazão constante. Neste caso um frasco de Mariotte fornecerá uma vazão constante da solução salina. Deve- 
se preparar a solução salina como descrito anteriormente e introduzido no frasco de Mariotte. 
Determinada a condutividade elétrica inicial. Daí o procedimento de medição da condutividade elétrica é 
exatamente o mesmo do já descrito. 


": 
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1.4.4. TRATAMENTO DOS DADOS OBTIDOS EM CAMPO 


Todos os dados devem ser organizados em planilha eletrônica Excel ou equivalente e trabalhados para a 
determinação da vazão. Fiquem atento com a condição de maré, pois se as medições levarem muito tempo 
podem haver variações nas condições de escoamento e consequentemente na vazão. 


Método do Flutuador: 


1. Calcule as três vezes a velocidade superficial; 

2. Faça a média das velocidades superficiais; 

3. Usando os coeficientes de correção da velocidade superficial com a velocidade média e determine 
o intervalo da velocidade média; 

4. Calcule a área da seção de escoamento; 

5. Determine a vazão. 


Método do Molinete: 


1. Determine a velocidade média em cada seção considerada; 
2. Calcule a área de influência de cada seção; 
3. Ache a vazão que passa por cada seção e a vazão total. 


Método da Diluição Química: 


Com os registros de concentração em função do tempo, volume de solução salina e concentração da 

Ci e su Rd E —  Ci-ch)Vi 

solução salina, usado a equação vista em sala, é possível determinar a vazão - Q O E a 
Jo te(b)=cnjdt 

1.4.5. COMPARAÇÃO DOS RESULTADOS 


Como serão usados métodos diferentes com variados procedimentos, fora a comparação direta dos 
resultados, também deve ser considerado a precisão entre as medidas, a reprodutibilidade entre os 
ensaios, as incertezas das medidas e aplicabilidades de cada método. 


Discutir as vantagens e desvantagens de cada método. 
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2.1. OBJETIVOS 


Esta aula tem por objetivo apresentar uma nova possibilidade para medição de vazão a partir de imagens: 
o uso de drones. 


Ao final da aula os alunos devem ter noção sobre os equipamentos e metodologia de medição aplicada 
durante a atividade desenvolvida. 


2.2. MEDIÇÃO DE VAZÃO E O USO DE DRONES 


Estudos envolvendo a utilização de Veículos Aéreos Não Tripulados (VANTSs) ou drones para medição de 
vazão vêm surgindo nos últimos anos como uma alternativa de medição remota em ambientes de difícil 
acesso e em condições de eventos hidro-meteorológicos adversos. As principais aplicações envolvem a 
medição de velocidades superficiais de rios por meio de técnicas baseadas em imagens, como, por 
exemplo, velocimetria de imagem de partículas em grande escala (LSPIV), velocimetria de rastreamento 
de partículas (PTV) e a velocimetria de imagens espaço-temporais (STIV) (LUKÓ, 2016). Estes métodos 
são baseados na análise de gravações de vídeo da superfície da água, para as quais a contribuição dos 
VANTSs tem se mostrado essencial (e.g TAURO ET AL., 2016; DETERT & WEITBRECHT, 2015; PERKS ET 
AL., 2016; SILVA ET AL, 2019). 


Para que seja possível captar e coletar imagens para posterior medição de vazão são imprescindíveis 
requisitos como: 


= Determinação de Pontos de Controle e Posicionamento: obtenção por GPS de três pontos de 
controle no terreno, constituídos por qualquer objeto na imagem cujas coordenadas X, Y, Z no terreno 
sejam conhecidas. Os pontos de controle devem ser bem distribuídos na imagem de modo a formar um 
triângulo e não uma linha reta. 


= Planejamento de Voo: esta etapa é fundamental para a adequada captação de imagens da área em 
estuda. A partir desta é possível estipular altitude de voo; número de linhas de voo necessárias; distância 
entre as linhas de voo; distância no terreno entre exposições; número total de exposições requeridas. 


Ainda nesta etapa é importante verificar a velocidade do vento. Ventos fortes podem causar vários 
problemas indesejados, como exemplo o desvio do drone da rota planejada, prejudicando a captura de 
imagens e posterior processamento das mesmas. 


= Execução do Voo e Processamento das Fotografias 


De modo geral os processamentos realizados com as fotografias correspondem a sequência: avaliação das 
imagens obtidas, alinhamento das fotos, construção da nuvem de pontos, construção da malha, construção 
do ortomosaico. Para esta etapa softwares devem ser utilizados. 


2.5. PRÁTICA EM CAMPO 


O método de medição de vazão proposto consiste em uma adaptação e aprimoramento do método do 
flutuador proposto por PALHARES et al. (2007). Neste caso deve ser estabelecido o tipo de flutuador 
usado na prática de medição. Nesta aula prática será usada serragem como flutuador. Recomenda-se a 
leituras dos trabalhos de avaliação de traçadores em estudos hidráulicos (PAULO, 2015; SILVA et at, 
2019). 


O local escolhido para aplicação deve ser um trecho retilíneo do rio e com informações pertinentes à 
seção, tais como medições da seção transversal e estação linimétrica. O local deve ainda oferecer base para 
os experimentos e alocação de pontos de controle (PC). 


Para esta aula prática serão utilizados os seguintes recursos e equipamentos: 
= Dispositivo móvel com GPS (determinação dos pontos de controle); 
= Molinete (medição de velocidade e posterior comparação); 


= Serragem; 
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n VANT DJI Mavic PRO. 
É necessário que o drone esteja com o GPS calibrado para o local de voo. 


As medições de vazão no rio são realizadas através levantamentos aerofotogramétricos com o drone (DJI 
Mavic PRO) na região a jusante da do ponto de lançamento da serragem. Os passos adotados para o 
processo do levantamento aerofotogramétrico, seguindo metodologia proposta por Silva et. al. 2019, são 
descritos a seguir: 


a) disposição do objeto de escala e dos pontos de controle dentro do campo de visão do local de 
medição do VANT; 
b) posicionamento do VANT sobre do rio, com a câmera disposta em um ângulo perpendicular ao 


fluxo do escoamento com altura suficiente para captura de imagens entre as duas margens do curso 
d'água; 


c) lançamento da serragem; 


d) Gravação de vídeo (câmera do drone) acompanhando o lançamento da serragem, de modo a 
registrar o avanço das partículas ao longo da superfície no rio; 


Feita a gravação, a vazão, posteriormente, deve ser estimada por meio da análise das velocidades das 
partículas em um processo de integração com a área da seção transversal do trecho do rio avaliado (ver 
item 2.6). 


2.6. PROCESSAMENTO DE DADOS 


O processamento de dados segue metodologia proposta por Silva et. al, 2019. Para avaliação da 
velocidade superficial das partículas, primeiramente deve ser realizada a conversão dos vídeos em 
imagens sequencias por meio da utilização do software livre Video To JPG Converter 
(https://www.dvdvideosoft.com/pt/products/dvd/Free-Video-to-JPG-Converter.htm) e após, a seleção 
das fotos de modo a se obter fotos com intervalos de 1 segundo. A avaliação das imagens pode ser 
realizada com o software Arcgis. 


A análise da velocidade superficial por meio de linhas permite demarcar a frente das plumas de traçadores 
com linhas, buscando aproximar-se de um traçado de perfil de velocidade superficial. Esta metodologia 
consiste em acompanhar, junto das regiões mais relevantes (centro, margem esquerda e margem direita), 
o movimento das frentes de aglomerados de serragem. Posteriormente, junto da referência de escala e o 
conhecimento do tempo decorrido foi possível determinar a velocidade superficial instantânea média 
entre duas imagens por meio da Equação 1. A vazão da seção, por fim, é determinada com a multiplicação 
da velocidade média da seção com a área da seção transversal e, quando necessário, com um coeficiente de 
ajuste (Equação 2). 


A TA (1) 


Onde: v é velocidade do traçador (m/s); d é distância percorrida pelo traçador, entre duas imagens, na 
escala da imagem (m); Lr é o comprimento real do objeto de escala (m); Lf é o comprimento do objeto na 
escala da imagem (m); e At é o tempo transcorrido entre as duas imagens analisadas. 


Q=Cx+vxA (2) 


Onde: € é o coeficiente de ajuste da velocidade média superficial em velocidade média da seção do rio; v é 
a velocidade do traçador (m/s); e A é a área da seção transversal analisada (m?). 
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3.1. OBJETIVOS INSTRUCIONAIS 
Ao final desta aula prática de infiltração da água no solo você deverá ser capaz de: 


e Enunciar quatro possíveis métodos aplicáveis a determinação da taxa de infiltração da água 
no solo; 


e Realizar a medição da taxa de infiltração em campo por um dos métodos testados na prática 
de campo; 


e Identificar as diferenças e aplicabilidades de cada método; 


e Comparar os resultados obtidos por diferentes métodos identificando suas aplicabilidades. 


3.2. INFILTRAÇÃO DA ÁGUA NA HIDROLOGIA 


Como comentado a visão de engenharia dada a nossa disciplina de Hidrologia e Climatologia foca o ciclo 
hidrológico do ponto de vista quantitativo. O conhecimento da taxa de infiltração da água no solo em 
hidrologia é fundamental para aplicar técnicas de conservação do solo, planejamento e dimensionamento 


de estruturas de infiltração e de drenagem, tal conhecimento aplicado ao planejamento urbano é 
extremamente útil no zoneamento municipal. 


Assim, para se dispor de valores de referência é fundamental o conhecimento dos métodos de medição e a 
espacialização dos dados como auxiliar à tomada de decisão. Em muitos casos o mapeamento dos 
diferentes tipos de solo associados a sua respectiva taxa de infiltração e as condições de saturação do solo 
são suficientes para tal aplicação. Portanto o conhecimento e da taxa de infiltração e seus interferentes é 
uma ferramenta extremamente útil na Engenharia Sanitária e Ambiental. Belém e em boa parte do Estado 
do Pará há uma peculiaridade que é a grande quantidade de sistemas de infiltração de efluentes, sendo a 
taxa de infiltração uma variável de interesse poucas vezes determinada. 


3.3. MÉTODOS DE MEDIÇÃO DE INFILTRAÇÃO DA ÁGUA NO SOLO 


Como discutido em sala de aula há diversos métodos para determinar a taxa de infiltração da água no solo, 
cada um com condições de operacionalização e ferramental específico, bem como diversas aplicabilidades 
com ênfase na área agrícola e zoneamento dos espaços urbanos. Como o foco da nossa aula é prático serão 
considerados os métodos descritos a seguir por se dispor de ferramental para sua execução. 


Os métodos de determinação da infiltração de água no solo são simples e capazes de representar 
adequadamente as condições do solo (PRUSKI et al., 1997). 


3.3.1. IFILTRÔMETRO DE ANÉIS CONCÊNTRICOS 


Para determinação da infiltração de água no solo o método mais conhecido é o infiltrômetro de anel (ou 
anéis concêntricos), onde a infiltração se processa apenas na vertical. Este consiste em dois anéis 
concêntricos (com 2 diferentes diâmetros), ambos de mesma altura, que são inseridos no solo até a 
profundidade de 0,15 m, com auxílio de marreta, na Figura 1 pode-se observar os componentes do 
infiltrômetro de anel, composto por 2 cilindros: 9500 x 250 mm e 9250 x 250 mm, espessura de 3 mm. 
Flutuador tipo boia com escala, suporte do flutuador, cruzeta em aço para cravação, nível tipo bolha, 
cronômetro digital portátil, marreta com pontas de Nylon e par de alças para extração. Sendo de suma 
importância que as bordas inferiores dos anéis sejam finas e apresentem corte em forma de bisel, para 
melhor introdução destes no solo. 
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Figura 1 - Componentes do infiltrômetro de anel. 


Fonte: https: //en.eijkelkamp.com/products /field-measurement-equipment/double-ring- 
infiltrometer.html 


Vale ressaltar que o ensaio de infiltração pode ser realizado na superfície do solo, como o que faremos, 
entretanto, pode acontecer em profundidades diferentes como mostrado na Figura 2. 


Figura 2 - Ensaio de infiltração com anéis concêntricos sendo realizado a um metro de profundidade da 
superfície do solo. 


Fonte: http://earthy-moony.blogspot.com/201 1/02 /measurement-of-infiltration-rate.html. 


No caso do dimensionamento de sumidouros como dispositivo de infiltração de efluentes deve-se garantir 
que a profundidade investigada seja representativa da profundidade onde se assentará a base de 
infiltração do sumidouro. A título de curiosidade, os dois recipientes de PVC da Figura 2 são frascos de 
Mariotte usados para manter o nível da água constante, isto é, mesma carga hidráulica. 
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3.3.2. ENSAIO DE INFILTRAÇÃO ABGE (1996) 


O método padronizado mais simples disponível para ensaios de infiltração é o recomendado pela 
Associação Brasileira de Geologia de Engenharia (ABGE) em 1996. 


O teste consiste em realizar um furo de sondagem com geometria cilíndrica no solo com o auxílio de uma 
cavadeira articulada. O furo precisa ficar com circunferência uniforme que pode ser 20 cm de diâmetro e 
atingir uma profundidade de 30 cm aproximadamente. O furo deve ser totalmente preenchido com água 
por 10 min garantindo que o solo ao redor esteja saturado, se necessário deve ser feita a reposição da 
água. 


Após os 10 min iniciais o ensaio pode começar. Neste momento efetua-se a contagem do tempo das 
medições de rebaixamento do nível d'água no furo, onde a mensuração do nível de água no furo de 
sondagem foi feita, ver Figura 3. Os intervalos, em minutos, pré-determinados e consecutivos sugeridos 
são: 0,5; 1; 1,5; 2; 2,5; 3; 3,5; 4; 5; 6; 7; 8; 10; 12; 14; 16; 18; 20 e 22 - período do ensaio geralmente 
previsto, procedendo-se as respectivas anotações em planilha (Apêndice 1), com o registro do decréscimo 
do nível d'água no furo, a cada intervalo de tempo informado. O ensaio pode ser encerrado quando o 
rebaixamento atingiu 20% da carga inicial ou 30 min. 


Figura 3 - Ilustração do ensaio de infiltração indicado por ABGE (1996) 


Para H de 0,20 a 0,30 em - ide 3a 5 


À equação a seguir é utilizada para a determinação da condutividade hidráulica de acordo com o método 
ABGE (1996): 


AH il 


Onde K (cm/s), é a condutividade hidráulica, AH (cm) é a variação do nível de carga hidráulica, At é a 
variação do tempo (s), H é a profundidade do poço (cm), r é o raio do furo do poço (cm) e i é adimensional 
variando de 3 a 5 (quando a geometria do furo é cilíndrica adota-se 5). 


3.3.3. MÉTODO OPEN END HOLE 


O método Open end Hole ou “Ensaio de Ponta Aberta” consiste na utilização de um furo de sondagem 
revestido com tubo de PVC, do topo até a sua base, para a obtenção do coeficiente de condutividade 
hidráulica vertical de uma determinada profundidade de solo (FIORI, 2010). Na Figura 4 apresenta-se as 
variáveis utilizadas para o ensaio. 
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Figura 4 - Ilustração do ensaio Open end Hole 
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A carga hidráulica (ho) é o comprimento do tubo, do topo até a base do tubo e é dado pela expressão hç = P 
+ L. O método consiste na avaliação da infiltração vertical carga hidráulica e/ou nível de água variável. 
Essa variação do nível de água (Ah) é dada em função do tempo, sendo assim h, representa o nível de água 
inicial e h; representa o nível de água final. Onde h; é dado pela relação: h; = ho + Ah. Considerando, ainda, 
o raio do tubo de PVC (R =D/2) e representa a condutividade hidráulica vertical: 


K, = 2,303 ( É ) | (2º)] 
ss a-At) [lh 


Onde Kv é a condutividade hidráulica vertical, R é raio, At é a variação do tempo, h, é a carga hidráulica 
inicial e h; é a carga hidráulica final. 


O ensaio é realizado em três pontos distintos espacializados de forma triangular (equilátera) distanciados 
em aproximadamente 1,5 m (Figura 5). Os furos de sondagem precisam ficar em diferentes profundidades 
0,5; 1,0 e 1,5 metros. Em locais onde a zona saturada é superficial, como Belém, atenção deve ser dada se 
fluirá água natural ao furo. Caso ocorra água, deve-se descartar o ponto e fazer o ensaio com os três furos, 
no entanto, com duas profundidades distintas (0,5 e 1,0 m), em caso de a um metro de profundidade 
também for observada a presença de água deve-se adotar três furos com a mesma profundidade (0,5 m). 
Em caso de ainda assim for constatada a presença de água com a 0,5 m então este ponto deve ser 
descartado e selecionado outro ponto para a realização do ensaio. Note que as profundidades 0,5; 1,0 e 1,5 
m não são obrigatórias e poderia ser adotado outro conjunto, por exemplo, 0,2; 0,4 e 0,6 m. 


Figura 5 - Disposição dos furos de sondagem para o ensaio Open end Hole 


Nível do solo 
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Um dos cuidados executivos para a realização do ensaio é que o tubo em PVC seja de diâmetro adequado 
para o furo de sondagem evitando caminho preferencial entre o a parede do tubo de PVC e a parede do 
furo de sondagem. 


3.4. PRÁTICA EM CAMPO 


Nossa aula prática acontecerá em campo fazendo medições reais, por diferentes métodos e com 
possibilidade de comparação dos métodos e resultados. 


3.4.1. SELEÇÃO DO LOCAL 


Para nossa aula de campo não há critérios específicos para a seleção do local de ensaio, pois trata-se 
apenas de uma aplicação prática. Todavia em aplicações na área agrícola deve ser usado algum critério 
para a divisão da área de interesse em subárea; em aplicações de zoneamento urbano deve ser 
identificados os tipos de solo e sua distribuição espacial, bem como evitar ensaios em locais de deposição 
de aterro; para o caso de infiltração de efluentes deve-se considerar a profundidade da zona saturada e a 
profundidade que se pretende posicionar o sumidouro. 


Para facilitar nossa atividade foi escolhida a lateral do prédio do Laboratório de Engenharia Sanitária e 
Ambiental (LAESA), no Setor Profissional, da Cidade Universitária Prof. José da Silveira Netto da 
Universidade Federal do Pará em Belém. Na Figura 6 é indicada a localização da lateral do prédio LAESA 
para a realização da aula prática. Na referida aparecem troco de tubo em amarelo que são as tampas de 
poços de monitoramento da água subterrânea, trata-se de estrutura de pesquisa e não deve ser usado 
como apoio ou para sentar! 


Figura 6 - Identificação da lateral do prédio da Faculdade de Engenharia Sanitária e Ambiental. 


3.4.2. ANTES DA MEDIÇÃO 

Separar os seguintes materiais antes de se deslocar para o local de medição: 
12) Protetor solar; 

13) Tenda ou guarda sol e cadeiras; 


14) Trenas de 3 e 50 metros; 


15) Medidor de nível d'água eletrônico; 

16) Planilhas de campo (Apêndice 1), prancheta, caneta e lapiseira; 
17) Caixa com todo o conjunto do infiltrômetro de duplo anel; 

18) Cavadeira articulada e cavadeira simples; 


19) Tubos de PVC de 100 mm, em 3 tamanhos; 
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20) Régua metálica. 


Não esqueça de estar vestido adequadamente de calça, botas, camisa confortável, óculos de proteção solar 
e chapéu/boné. Também fiquem atentos para fazer registros fotográficos em todas as etapas. 


3.4.3. ATIVIDADES EM CAMPO 


Em grupos de 5 alunos nesta atividade de campo serão testados os três métodos de medição da taxa de 
infiltração e consequente condutividade hidráulica do solo. 


Inicialmente um grupo fará medição usando o primeiro método, concomitante em seguida outro dois 
grupos farão a medição usando o segundo e terceiro métodos. Caso haja tempo disponível todos os grupos 
farão a medição dos seus próprios dados, caso contrário os dados medidos por cada grupo deveram ser 
compartilhados com os demais grupos para a elaboração do seu respectivo relatório. 


Os ensaios podem ser iniciados ao mesmo tempo em pontos distintos guardando alguma distância para 
não atrapalhar a outra equipe sendo: 


1) Método ABGE (1996): aqui será considerado o Ensaio em Poço, sem revestimento e acima do nível 
d'água. Usando a cavadeira articulada faça um furo de aproximadamente 30 cm de profundidade. Meça o 
diâmetro aproximado do furo e a profundidade que ficou. Sature o furo com água até a superfície por 10 
mim. Passados os 10 min complete o furo com água e comece a fazer os registros da variação do nível da 
água com auxílio de uma régua que deve ficar dentro do furo, considerando os intervalos indicados na 
planilha (Apêndice 1). Segundo ABGE (1996) o ensaio pode ser encerrado quando o rebaixamento atingiu 
20% da carga inicial ou 30 min. 


2) Infiltrômetro de duplo anel: abra a caixa de transporte do equipamento e retire os anéis, cuidado com o 
peso da caixa e dos anéis. Posicione os anéis metálicos na área de interesse, lembrando de colocar a lateral 
mais fina de cada anel para baixo. Ajuste a posição dos anéis com a ferramenta de ajuste no formato de 
uma cruzeta. Use marreta da caixa para fazer a cravação de aproximadamente 10 cm dos anéis no solo. 
Consulte se estão nivelados com a régua de nível, faça eventuais ajustes. Coloque inicialmente água no anel 
externo até aproximadamente formar uma lâmina de 10 cm. Posicione o suporte e a régua do cilindro 
central. Acrescente água até formar uma lâmina de 10 cm. Marque o nível inicial e comece a contagem do 
tempo. Faça a anotação dos registros de altura da lâmina d'água nos intervalos de tempo indicados na 
planilha de campo (Apêndice 1), caso o nível rebaixe mais que dois centímetros faça a reposição da água 
nos anéis. Siga a orientação das colunas da planilha de campo. Depois de duas leituras de rebaixamento é 
possível calcular demais colunas, sendo a terceira coluna a diferença entre as medidas, a quarta coluna os 
infiltrados acumulados e a quinta coluna a velocidade de infiltração, atenção as velocidades de infiltração 
são em mm/h, então será necessário multiplicar o tempo medido em minutos por 60. O ensaio continua 
até que sejam atingidas três velocidades de infiltração iguais de forma consecutiva. 


3) Método Open end hole: o ensaio será realizado em 3 profundidades distintas, sendo elas: 0,3; 0,6 e 0,9 m. 
Faça 3 furos com cavadeira articulada ou trado helicoidal nas profundidades indicadas, seguindo uma 
forma triangular (equilátera), com lado de 1,5 metros. Encaixe os tubos de PVC em cada um dos 
respectivos furos, lembre-se que o furo precisa estar ajustado ao tubo para não escapar água entre o furo e 
o tubo de revestimento. Faça leituras regulares do nível d'água a cada 5 min por um intervalo de ao menos 
40 min. Preencha as informações da planilha de campo (Apêndice 1). 


3.4.4. TRATAMENTO DOS DADOS OBTIDOS EM CAMPO 


Todos os dados devem ser organizados em planilha eletrônica Excel ou equivalente e trabalhados para o 
cálculo da taxa de infiltração e da condutividade hidráulica. 


Método ABGE (1996): 
6. Calcule as taxas de infiltração; 


VÊ Faça gráficos da variação do nível d'água e da taxa de infiltração em função do tempo; 
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8. Com a equação fornecida calcule a condutividade hidráulica da zona não saturada do solo. 
Infiltrômetro de duplo anel: 
4. Calcule as taxas de infiltração; 


5. Faça gráficos da variação do nível d'água e da taxa de infiltração em função do tempo. 


Método Open end hole: 


1. Calcule as taxas de infiltração; 
2 Faça gráficos da variação do nível d'água e da taxa de infiltração em função do tempo; 
3. Com a equação fornecida calcule a condutividade hidráulica vertical da zona não saturada do solo. 


3.4.5. COMPARAÇÃO DOS RESULTADOS 


Como serão usados métodos diferentes com variados procedimentos, fora a comparação direta dos 
resultados, também deve ser considerado a precisão entre as medidas, a reprodutibilidade entre os 
ensaios, as incertezas das medidas e aplicabilidades de cada método. 


Discutir as vantagens e desvantagens de cada método. 


Apresente uma aplicação prática dos resultados para o dimensionamento de uma trincheira de infiltração 
de água da chuva da cobertura de uma casa ou prédio e para um sumidouro que infiltrará efluentes 
tratados de uma casa ou qualquer outro empreendimento. 
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4.1. OBJETIVOS INSTRUCIONAIS 

Ao final desta aula prática de Estatística Aplicada à Hidrologia você deverá ser capaz de: 
* Acessar ao sítio da Agência Nacional de Águas para localizar e copiar dados hidrológicos; 
e Enunciar quatro distribuições probabilísticas aplicadas frequentemente em Hidrologia; 


e Testar a aderência das diferentes distribuições probabilísticas a dados medidos usando 
planilha eletrônica Excel; 


e Identificar que tipo de programa específico pode ser utilizado na análise de frequência. 


4.2. ESTATÍSTICA APLICADA A HIDROLOGIA 


Séries de variáveis hidrológicas como precipitações, vazões, nível d'água, evaporação e outras, quando 
observadas ao longo do tempo, apresentam variações sazonais. Estas variações não são regulares, isto é, 
ocorrem anos com valores máximos bem acentuados e outros anos com valores reduzidos. A observação 
de séries longas de dados hidrológicos revelará a ocorrência de extremos (máximos e mínimos) e 
diferentes sequências de valores, que caracterizam as variáveis hidrológicas como aleatórias. Variáveis 
hidrológicas estarão sempre associadas, portanto a uma probabilidade de ocorrência. Consequentemente 
as obras hidráulicas devem sempre ser dimensionadas para um determinado “risco” de falha. 


Neste caso o objetivo da estatística é o de extrair informações significativas de uma dada massa de dados. 
As técnicas utilizadas em estatística aplicadas à hidrologia permitem avaliar a probabilidade de ocorrência 
de um fenômeno hidrológico com determinada magnitude. 


4.3. EXTREMOS MÁXIMOS: PREVISÃO DE ENCHENTES 


Quando se trabalha com valores máximos de uma dada variável hidrológica estão sendo considerados os 
extremos máximos. No caso da previsão de enchentes são trabalhadas as séries de dados de vazão, 
selecionando-se apenas as vazões máximas de cada ano em um determinado local. Neste caso é obtida a 
série de vazões máximas deste dado local e é possível realizar análises estatísticas relacionando vazão 
com sua probabilidade associada. As séries de vazões disponíveis na maior parte dos locais (postos 
fluviométricos) são relativamente curtas, não superando algumas dezenas de anos. Por outro lado, muitas 
vezes é necessário a estimativa de uma vazão para uma condição futura desconhecida, que pode estar 
associada uma centena de anos ou até século. 


4.4. LEVANTAMENTO DE DADOS 


A elaboração de projetos em Hidrologia e Recursos Hídricos, por vezes, requer o levantamento de 
informações. Tais informações podem ser obtidas em levantamentos em campo, que são demorados e 
apresentam um custo, mas de grande valia conforme a necessidade de detalhamento do projeto. 


Outra fonte de informações com dados é Sistema Nacional de Informações sobre Recursos Hídricos 
(SNIRH), gerenciado pela Agência Nacional de Águas (ANA), é parte integrante dos instrumentos de gestão 
previsto na Política Nacional de Recursos Hídricos, instituída pela Lei nº 9.433, de 08 de janeiro de 1997, 
conhecida como Lei das Águas. Trata-se de um amplo sistema de coleta, tratamento, armazenamento e 
recuperação de informações sobre recursos hídricos, bem como fatores intervenientes para sua gestão. 


No presente material será considerado o uso de dados disponíveis no Portal HidroWeb que é uma 
ferramenta integrante do SNIRH e oferece o acesso ao banco de dados que contém todas as informações 
coletadas pela Rede Hidrometeorológica Nacional (RHN), reunindo dados de níveis fluviais, vazões, 
chuvas, climatologia, qualidade da água e sedimentos. 
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4.5. PRÁTICA EM LABORATÓRIO DE INFORMÁTICA 


Como dito os dados usados nesta prática serão obtidos do Sistema HIDROWEB gerenciado pela ANA. Para 
tanto acesse ao sítio: 


- http: //www.snirh.gov.br /hidroweb/apresentacao 
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- Nesta nova página você poderá optar por inserir tipo de estação (pluviométrica ou fluviométrica ou 
ambas), o código da estação desejada, ou o nome da estação, ou o nome do rio (precisa selecionar a Bacia 
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“Consultar”. Como resultados da consulta vão aparecer para todas as estações contidas na bacia 
selecionada. 
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- No retorno da busca, parte de baixo da tela, é possível visualizar uma lista de estações, com o código, 
nome e tipo de estação. 
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Por continuidade atribui-se a denominação de Bacia 9 à área de drenagem de qualquer bacia hidrográfica 
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monitoramento, a grande maioria dos campos é constituída por zeros (000). 
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- Daí é só você localizar a estação do seu interesse, no campo “Tipo de arquivo” marcar um dos três 
formatos e clicar em “BaixarArquivo” que o download será iniciado. 
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- Seu arquivo de dados está pronto para ser trabalhado. 
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4.5.1. PREPARAÇÃO DOS DADOS 


Nesta disciplina nos interessamos em separar as maiores vazões de cada ano. Desta forma você precisará 
separar as maiores vazões de cada ano. 


Como exemplo iremos usar as vazões máximas anuais do rio Guaporé, no posto fluviométrico (86560000) 
Linha Colombo (RS), que já foram separados e organizados e são apresentados na Tabela 1, a seguir. 


Tabela 1 - Vazões máximas anuais do posto fluviométrico (86560000) Linha Colombo (RS) 


Vazão Máxima 


Ano Ano Vazão Máxima (mº/s) Ano 


(m2/s) 


1940 
1941 
1942 
1943 
1944 
1945 
1946 
1947 
1948 
1949 
1950 
1951 
1952 
1953 
1954 
1955 
1956 
1957 598 1975 480 1993 


Vazão Máxima 


O) 


Falha 
673 
760 
780 
653 
537 
945 

1650 

1165 
888 
728 
809 
945 

1380 

Falha 

Falha 

Falha 
639 


Vale ressaltar que todos os anos onde ocorram falhas foram desconsiderados. Todavia poderia ser usado 
algum critério de preenchimento de falha, porém não é o objetivo da presente aula. 


o 
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4.5.2. DETERMINAÇÃO DE VALORES MÁXIMOS POR DIFERENTES DISTRIBUIÇÕES PROBABILÍSTICAS 


Todos os projetos de engenharia são planejados para o futuro, não havendo certeza absoluta das exatas 
condições de trabalho da obra. O engenheiro que atua com hidrologia e recursos hídricos não estará 
absolutamente certo da vazão ou precipitação que afetará o projeto. Contudo, deve estar consciente de que 
um erro acentuado de previsão das quantidades hidrológicas poderá causar efeitos destruidores 
indesejáveis, que podem inviabilizar economicamente todo o projeto. 


Uma vez que o comportamento exato das vazões em anos futuros não pode ser absolutamente previsto, 
procura-se introduzir leis de probabilidade de modo a estabelecer as prováveis variações para permitir 
que o plano seja completado com base em um risco calculado. Recorre-se, pois, à análise estatística com o 
propósito de utilizar os eventos de descargas observadas (série histórica de vazões) num dado período, 
como meio de se efetuar a projeção para um período de tempo maior. 


Na previsão de enchentes, ou seja, na determinação da magnitude das vazões de pico das cheias (que são 
as vazões críticas ou de projeto), recorre-se ao uso de modelos de probabilidade, a partir de um enfoque 
estatístico que consiste em definir a relação entre as descargas máximas e as correspondentes frequências 
de ocorrência, apoiando-se no estudo de uma série de dados observados. A suposição básica é que as 
cheias verificadas durante um determinado período possam ocorrer em um período futuro de 
características hidrológicas similares, isto é, com uma expectativa de repetição. 


As funções matemáticas de distribuição de probabilidade mais utilizadas na análise de frequência das 
vazões de enchente são: 


1) distribuição gama, também conhecida como distribuição Pearson tipo III; 


2) transformação logarítmica da distribuição gama, também conhecida como distribuição log-Pearson tipo 
HI; 


3) distribuições exponenciais, também conhecidas como distribuições de valores extremos ou 
distribuições de Fisher-Tippett, que são de três tipos, sendo a mais frequente do tipo I, duplo exponencial, 
conhecida como distribuição Gumbel; 


4) distribuição gaussiana (distribuição normal de probabilidade); 


5) transformação logarítmica da distribuição normal, também conhecida como distribuição log-normal ou 
distribuição de Galton. 


Ven Te Chow conseguiu comprovar que a maioria das distribuições de probabilidade usadas em 
Hidrologia por uma única equação geral, onde a variável de interesse (vazão, chuva, etc.) é expressa em 
função da média, do desvio padrão e do fator de frequência Kr;, conforme: 


XTr =D oo gu DS 
Onde: 
Xry = Variável de interesse (vazão, chuva, etc.) para o período de retorno Tr; 
X = média aritmética do conjunto de dados; 
Kry = fator de frequência associado a distribuição probabilística de interesse; 


Sxy= desvio padrão do conjunto de dados. 


Como visto na aula teórica em sala de aula, a média aritmética, o desvio padrão e coeficiente de assimetria 
são determinados diretamente da série histórica dos dados que pode ser determinado diretamente com 
calculadora, planilha eletrônica (função do Excel: “=média()”, “=desvpad.a()” e “=distorção()”) ou software 
específico. No caso das distribuições que trabalham com o logaritmo (log-normal, log-Pearson, etc.) é 
necessário antes de calcular a média aritmética, o desvio padrão e o coeficiente de assimetria determinar o 
logaritmo (pode ser logaritmo na base 10, log, ou logaritmo natural, In, usaremos a segunda opção) de 


cada valor da variável de interesse (vazão, precipitação, etc.). 
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O valor do fator de frequência pode ser obtido por tabela, como indicado na apostila e feito em sala de 
aula. Também pode ser usado o Excel, neste caso, para as distribuições Normal e log-Normal, o Krr, é dado 
por: “= INV.NORMP.N(1-(1/Tr); Média; Desvio Padrão)”. 


Para a distribuição Gumbel segue a mesma ideia usada em sala de aula, isto é, calcula-se a variável 
reduzida de Gumbel (yrr) para cada período de retorno, em seguida determina-se a média e desvio padrão 


dos valores de “yr” e finalmente aplica-se a equação: K YTr-Ymédio 


Trgumbel=5 
desvio padrão 

No caso das distribuições Pearson III e log-Pearson III, que são derivados da função Gama, usa-se a função 

do Excel “= INV.GAMA(1-(1/Tr); «; B)”. Todavia precisa feito um ajuste na equação proposta pode Ven Te 

Chow. 


Inicialmente, sendo o coeficiente de assimetria positivo (gx>0), os parâmetros “a e B” da distribuição Gama 
são dados por: 


4 
4A=— 
2 
Ix 

Onde: gx - coeficiente de assimetria da variável de interesse (vazão, chuva, etc.). 
Sx 


à qa 


Onde: Sx - desvio padrão da variável de interesse (vazão, chuva, etc.). 


Para o Excel a equação da distribuição Pearson III fica no formato: 


2º Sy 


E 1 
Xrr = X — + INV.GAMA (1 <a; 6) 


Caso o coeficiente de assimetria positivo (gx<0), os parâmetros “x e B” da distribuição Gama são dados 
por: 


Neste caso, para o Excel, a equação para a distribuição Pearson III fica no formato: 


5x 
9 


E il 
Xrr = X — — INV.GAMA (= 1a; 6) 


Para a distribuição log-Pearson III o procedimento é parecido com o anterior, lembrando que para todas se 
as distribuições “log” inicialmente os valores da variável de interesse devem ser logaritmados na forma Y 

= In(X), os coeficientes parâmetros a e 8 precisarão ser novamente calculados, considerando a média 
aritmética, desvio padrão e coeficiente de assimetria dos In(X). 
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Para o Excel, caso gy>0, a equação da distribuição log-Pearson III fica no formato: 


Ze Soy 


E 1 
Xry = exp vá — ZE + INV. GAMA (1 Tr a; 3) 


No caso de “gy< 0” a equação da distribuição log-Pearson III fica no formato: 


Sy 


an; 1 
Xry = exp Dá — — INV. GAMA (=: ; )) 


Iy 


4.5.3. POSIÇÃO DE PLOTAGEM PARA ANÁLISE DE FREQUÊNCIA 


Em hidrologia existem muitas vezes eventos extremos, tais como secas ou cheias. O valor de um 
acontecimento extremo é inversamente proporcional à sua frequência de ocorrência, isto é, um 
acontecimento extremo ocorre com menos frequência do que um evento moderado. Neste caso aplica-se a 
análise frequência em hidrologia para relacionar a magnitude dos valores extremos com a sua frequência 
de ocorrência, daí é possível comparar tais valores extremos com os determinados por distribuições 
probabilísticas e em seguida fazer previsões para caso mais extremos. Os resultados desta análise podem 
ser usados em vários problemas de engenharia, tais como, dimensionamento de barragens, pontes, 
estruturas de controlo de cheias, etc. 


No caso de eventos máximos, que o objetivo deste material, a estimação da probabilidade empírica de 
excedência, associada a um certo ponto, é geralmente referida como a determinação da posição de 
plotagem, a qual pode ser expressa como uma fração entre O e 1, ou como uma porcentagem entre O e 100. 
No caso da população, a determinação da posição de plotagem é meramente um problema de 
determinação da “fração” da população cujos valores são maiores ou iguais ao valor em questão. Assim, 
para uma variável aleatória de máximos, o menor valor da população terá uma posição de plotagem igual a 
1 (um) e o maior valor terá uma posição de plotagem igual a O (zero). A definição das posições de 
plotagem para dados amostrais não é tão direta como no caso populacional, uma vez que nunca haverá 
certeza de que a amostra contém o maior e/ou o menor valor da população. Assim, para dados amostrais, 
as posições de plotagem 0 e/ou 1 devem ser evitadas, à exceção dos casos em que existirem informações 
definitivas acerca dos limites populacionais (Naghettini e Pinto, 2007). 


Em síntese, para o caso de máximos, uma fórmula para a estimação da posição de plotagem deve 
especificar a frequência com que um, entre n valores ordenados de modo decrescente, será igualado ou 
superado. 


As posições de plotagem mais usuais em Hidrologia são dadas na Tabela 2. Lembrando que p = 1/Tr. 


Tabela 2 - Equações para estimar a posição de plotagem 


Nome Equação 
Weibull E 
Edi E +1 
, m-— 0,3175 
Mediana p=—————>— 
n+ 0,365 
— 0,35 
Hosking p= ia 
n 
m-— 0,375 
Blom p=——— 
n + 0,25 
m-—0,4 
Cunnane p= 
n+ 0,2 
: m — 0,44 
Gringorten =—D——— 
n + 0,12 
m-—0,5 
Hazen = : 
á n 
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Obs: quando os dados são ordenados em ordem decrescente o número de ordem “m” varia de 1, maior 


Aus 


valor, até “n”, que será o último valor. 


Não existe uma posição de plotagem melhor ou pior, em geral observa-se a melhor disposição dos dados. A 
posição de Plotagem de Weibull é indicada pelo US Water Resources Council. Já a publicação britânica 
Flood Estimation Handbook, mencionada por Bayliss e Reed (2001), sugere a utilização da equação de 
Gringorten. 


Aqui no nosso curso de Climatologia e Hidrologia da UFPA será usada a posição Weibull como referência. 


4.5.4. TESTE DE ADERÊNCIA 


Por hipótese assumimos que a distribuição populacional da variável de interesse (vazões, chuva, etc.) seja 
descrita por determinado modelo de distribuição de probabilidade (normal, log-normal, Gumbel, etc.). 
Testaremos se a hipótese verificando a aderência dos dados da amostra (série de vazões, chuva, etc.) ao 
modelo de distribuição de probabilidade. A hipótese de nulidade Ho é, sempre, que a distribuição teórica 
representa adequadamente a distribuição observada. 


Neste momento vamos usar o teste Não Paramétrico de Hipótese de Kolmogorov-Smirnov para conferir a 
aderência das distribuições probabilística aos dados. 


O Teste de Aderência de Kolmogorov-Smirnov permite verificar a significância das diferenças entre as 
frequências observadas e as frequências teóricas obtidas a partir de uma distribuição probabilística. 
Sendo aplicável para distribuições probabilísticas contínuas. 


Definidas as funções probabilísticas acumuladas empírica S(x) (observadas, isto é, no nosso caso, dados 
copiados do Hidroweb associados a uma posição de plotagem) e teórica F(x) (hipotética, distribuições 
probabilísticas testadas no Excel), a estatística de teste é o máximo módulo da diferença das 
probabilidades acumuladas destas funções: D = máximo|F -S], o valor de “D” é comparado com o valor de 
“Do” obtido da Tabela 3, sendo “N” o número de dados e “o” a significância. Se “D > Do” a hipótese da 
distribuição teórica representar a distribuição empírica é rejeitada ao nível de significância a. Se “Da > D” 
a hipótese da distribuição teórica representar a distribuição empírica não pode ser rejeitada ao nível de 


significância q. 


Tabela 3 - Valores críticos da estatística Dya do teste de aderência Kolmogorov-Smirnov 


10 0,369 | 0409 | 0,457 | 0,489 26 0,233 0,259 0,290 0,311 
mi 0,352 | 0,391 | 0,437 | 0,468 E 0,229 0,254 0,284 0,305 
12 0,338 | 0,375 | 0,419 | 0,449 28 0,225 0,250 0,279 0,300 
13 EEE ER 29 0,221 0,246 0,275 0,295 
14 0,314 | 0349 | 0,390 | 0,418 30 0,218 0,242 0,270 0,290 
15 0,304 | 0,338 | 0,377 | 0,404 31 0,214 0,238 0,266 0,285 
16 0,295 | 0327 | 0,366 | 0,392 32 0,211 0,234 0,262 0,281 
7 2660 DOG Bossa Dos 33 0,208 0,231 0,258 0,277 
18 0,279 | 0,309 | 0,346 | 0,371 34 0,205 0,227 0,254 0,273 
19 EATON DEE DUE E 0,202 0,224 0,251 0,269 
20 0,265 | 0,294 | 0,329 | 0,352 36 0,199 0,221 0,247 0,265 
21 2500 [0207] [OS DIGA 37 0,196 0,218 0,244 Ea 
22 0,253 | 0281 | 0,314 | 0,337 38 0,194 0,215 0,241 0,258 
23 AA Dodo soa E 39 0,191 0,213 0,238 0,255 
24 0,242 | 0,269 | 0,301 | 0,323 40 0,189 0,210 0,235 0,252 
LE 0,238 | 0,264 | 0,295 | 0,317 >40 | 1,22/raiz(N) | 1,36/raiz(N) | 1,52/raiz(N) EA 
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DN) 


No nosso curso iremos adotar um nível de significância “a” de 5% ou 0,05. Tal valor é muito empregado na 
literatura sobre Estatística aplicada a Hidrologia. 


O teste Kolmogorov-Smirnov pode ser implementado em planilha eletrônica Excel com os seguintes 
passos: 


1) Ordene os dados da sua variável de interesse do maior para o menor (considerando que estamos 
trabalhando com valores máximos), calcule o também os “In” da variável de interesse, determinando a 
média aritmética, o desvio padrão e o coeficiente de assimentria; 


2) Atribua um número de ordem “m” variando de 1 associado ao maior valor da variável de interesse 


“un 


até “n” que está associado ao menor dos valores; 


3) Usando a posição de plotagem de Weibull calcule o período de retorno associado aos dados. Em 
seguida calcule a probabilidade excedência fazendo “pexc = 1/Tr”. Depois calcule a probabilidade de 
inexistência “Pinex = 1-Pexc”; 


4) Como o teste de Kolmogorov-Smirnov pede que seja calculada diferença entre a probabilidade de 
inexistência da variável de interesse com probabilidade de inexistência de cada distribuição probabilística 
considerada: 


a. Distribuição Normal: use a “=DIST.NORM.N (Variável de interesse; Média; Desvio Padrão; 
VERDADEIRO); 
b. Distribuição log-Normal: use a “=DIST.LOGNORM.N (In(Variável de interesse); Média In(x); Desvio 


Padrão In(x); VERDADEIRO); 


Cc. Distribuição Gumbel: calcule os parâmetros “u = média - 0,577*c”, “a = 0,7797*Sx” e “=EXP(- 
EXP((2-Variável de interesse)/0))”; 
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d. Distribuição Pearson II: calcule os parâmetros “a = (2/gx)?”, “B = Sx/raiz(o)”, 'y =X —-Sy: Va” e 
“=DIST.GAMA(Variável de interesse -y; o; 8; VERDADEIRO)”; 


e. Distribuição log-Pearson III: calcule os parâmetros “a = (2/gy)”, “B = Sv/raiz(o)”, y=Y — Sp: Va” 
e “=DIST.GAMA(In(Variável de interesse) -y; «; 8; VERDADEIRO)”. 


5) A diferença das probabilidades acumuladas para cada distribuição probabilística é calculada 
como: D = máximo|pinex- Pinex(distribuição de interesse)|; 


6) Todas as distribuições probabilísticas que satisfizerem a condição “Da > D” suportam a hipótese 
da distribuição teórica não pode ser rejeitada ao nível de significância considerado. 


4.5.5. SUA ATIVIDADE PRÁTICA 
1) Trabalhando em grupos de até 5 alunos; 


2) Escolha uma estação, pode ser Pluviométrica ou Fluviométrica, contida no Sistema HIDROWEB 
gerenciado pela ANA que contenha pelo menos 20 anos de registros que podem ser contínuos ou 
descontínuos; 


3) Identifique o nome e o número da estação considerada; 

4) Separe os maiores valores da variável de interesse de cada ano; 

5) Ordene os dados do maior para o menor valor; 

6) Atribua um número de ordem, calcule a frequência (posição de plotagem de Weibull) e o período 
retorno; 

7) Calcule o logaritmo neperiano de cada valor da variável de interesse; 

8) Determine a média aritmética, o desvio padrão e coeficiente de assimetria da variável de interesse 


e dos respectivos logaritmos neperianos; 


9) Faça o cálculo comparativo entre os valores da variável de interesse com os fornecidos pelas 
distribuições probabilísticas: Normal, log-Normal, Gumbel, Pearson II e log-Pearson III. Considerar os 
mesmos períodos de retorno já calculados; 
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10) Faça um Gráfico comparativo entre os valores da variável de interesse com os fornecidos pelas 
distribuições probabilísticas em função dos respectivos períodos de retorno; 


11) Por comparação gráfica indique qual distribuição probabilística apresenta a melhor aderência aos 
dados da variável de interesse; 


12) Bônus para quem aceitar o desafio: aplique o teste de Kolmogorov-Smirnov e indique quais das 
distribuições probabilísticas testadas podem ser aceitas. 


4.5.6. USANDO O PROGRAMA SISTEMA ESPECIALISTA EM ANÁLISE DE FREQUÊNCIA 


Dentre as ferramentas disponíveis para a análise de frequência local, os sistemas especialistas 
computacionais, que emulam os princípios de raciocínio de um especialista humano ao selecionar uma 
distribuição de probabilidades, têm-se mostrado muito úteis. Um exemplo desse tipo de sistema é o SEAF 
(Sistema Especialista para Análise de Frequência local de eventos máximos anuais). Vamos testar o 
programa SEAF (Sistema Especialista em Análise de Frequência) 1.0, desenvolvido pelo Departamentos de 
Engenharia Sanitária e Ambiental e Engenharia Hidráulica e Recursos Hídricos da Universidade Federal de 
Minas Gerais, disponível na URL http://www.ehr.ufmg.br/downloads. 


H - Exemplo 01 - Excel GiovamiPene E m —- o 


Arquivo Págins Inicial Inserir Layout da Página Fórmulas Dedos Revisão Exibir Ayuda ACROBAT 9 Diga-me o que você deseja fazer 2 Compartilhar 


) 3 | Beinserir - - 
22 Quebrar Texto Automaticamente | [Geral z Dj | Ef Inserir z by o) 
Fl d E Excluir » - 
Formatação Formatar como Estilasde | = Classificar Localizar e 
Condicional- Tabela  Célular | ElFormataro | É criirar- Selecionar” 


Estilos Células Edição 


[El Mesclar e Centralizar = E. 9% 00 “4 aê 


Alinhamento y Número 


exemplo 01 


EE É Digite aqui para pesquisar 


Basicamente os dados da variável de interesse só precisam ser colocados no formato que é lido pelo SEAF. 
Os dados devem estar em uma única coluna e salvos em .txt, tal operação pode ser facilmente realizada no 
Excel. Como mostrados nas Figuras abaixo. 
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D exemplo 01 - Bloco de Notes 
Arquivo Editar Formatar Exibir Ajuda 


Linha 1, Coluns 1 100% Windows (CRLF) uT-s 


£ Digite aqui para pesquisar à ezasS ag O O Rd] iam q 


Em seguida o programa SEAF pode ser executado. No meu caso a versão do SEAF disponível não roda nas 
versões atuais do Windows, daí foi usada uma Máquina Virtual, Oracle VM VirtualBox, para rodar o 
Windows XP, conforme Figura a seguir. 


4 Oracle VM VirtualBox Gerenciador 
Arquivo F) Máquina Ajuda (H) 


GU! revamena 2a o. 


Novo Configurações Descarta  Eridar (1) 
E) Geral 


Nome: win x26s 
Sstema Cpereconal: viindans X? (32-it) 


EX) sistema 


Memória Pnnaasl: 192M0 

Ordemde Boot: Disquete, Óptico, Disco Rigido 
Aceleração: VT=x/AMD-N, Paginação Arinhada 
E) visplay 

Memória de Vídeo: 1546 

Controladora Gráfica: VBOXVGA 

Sorvgor de Desktop Remoto: Desabátado 
Gravação: Desabltado 

E) Armazenamento 


Controladora: IDE 

IDE Prméno Master: Win XPES.vci (Normal, 10,00 GE) 

IDE Secundáro Master: [Disco Óptco] vEcaGuestáccibans so (56, 76 MB) 
jo Audio 
Ortver do Hoggedero:  vindoas Drectsound 
Controladora: ICH ACS7 
E rede 
Adaptador 1: Intel PRO/1009 7 Server (NAT) 
O us 


Controladora USE: oscI 
Filtros de Dispestvo: O (0 atvos) 


Pastas Compartilhadas 
Pastas Compartihadas: 1 
E) Descrição 

Menus 


ES & Digite aqui para pesquisar ones. Yy o o dê Go) Rg [mt 


Agora é possível observar SEAF no Windows XP. 
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nm 


Meus documentos 


| 
| 


Meu computador 


SEAF iniciado. 


a 


Meus documentos 


Meu computador 


“à SEAF - SISTEMA ES 


Projeto Ferramentas Ajuda 
D|B/z| E 


SEAF - PROGRAMA PARA ANÁLISE DE FREQUÊNCIA DE EVENTOS 
HIDROLÓGICOS MÁXIMOS ANUAIS 


Versão 1.0 


Autor: Márcio de Oliveira Cândido 
Ena. Hidrólogo da CPRM 
E-mail: marciobh.cprm.gov.br 


wordPad 


Orientador: Mauro da Cunha Naghettini 
Professor da EEUFMG 


Fi-Ajuda F10-Opções 


———+ 


“5 Iniciar 


Projeto sendo criado. 
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El 


“ se” | Ocidental ” N Z 5 2 


Re 2 SA BG De BS Gsi io: 
“Eno Hictólooo da CPRM 


SPECIALISTA E ) 
ALISTA AP 


Ferramentas Ajuda 


CI 


IGRAMA PARA ANÁLISE DE FREQUÊNCIA DE EVENTOS 
HIDROLÓGICOS MÁXIMOS ANUAIS 


Versão 1.0 


dutor: Márcio de Oliveira Cândido 
Eng. Hidrólogo da CPRM 
E-mail: marcio&bh. cprm.gov.br 


Orientador: Mauro da Cunha Naghettini 
pepr=çi =” Professor da EEUFMG 


“3 Iniciar 
Fi-ajuda FilDO-Opções 


< 
Para obter ajuda, pressione Fi 


“s Iniciar 


La 


Meus documentos 


pl 


Pet computado! 


EAF - Criar projeto 


Título: [âula de Hidrologia Estatística 


Local: Belém 
Comentário: Arquivo do Professor 


Arquivo de dados: [Z2:Aseaftexemplo OM.txt 


&S Criar | [4.4] Cancelar | 


Fi-ájuda Fi0-Opções 


Iniciar 
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Depois dos dados lidos é só seguir pressionando a seta de execução que os resultados são apresentados. 


Ao final um arquivo “txt” com os resultados é apresentado. O descritivo do arquivo “txt” é apresentado a 
seguir: 


SEAF - SISTEMA ESPECIALISTA EM ANÁLISE DE FREQUÊNCIA 
TÍTULO: Aula de Hidrologia Estatística 

LOCAL: Belém 

COMENTÁRIO: Arquivo do Professor 

ARQUIVO DE DADOS: Z:Aseaflexemplo 0L.txt 


ESTATÍSTICAS DESCRITIVAS: 
**Tamanho da amostra: 47 
Máximo: 1785.0000 

Mínimo: 246.0000 

Média: 746.5532 

Desvio padrão: 350.7593 
Assimetria: 0.9869 


**Momentos-L: 
1: 746.5532 
2: 193.2720 
t3: 0.1807 

t4: 0.1535 


**Logaritmos: 
Máximo: 7.4872 
Mínimo: 5.5053 
Média: 6.5091 

Desvio padrão: 0.4730 
Assimetria: -0.1354 


TESTES NÃO PARAMÉTRICOS: 
**Teste de Kendall: 

Coef. Corr. Kendall: 0.1198 
Estatistica Z: 1.1650 

Valor limiar de | Z |: 1.9600 
Alfa: 5% 


Z<|Z | => dados independentes 


**Teste de Mann-Kendall: 
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Estatistica U: 1.6146 
Valor limiar de |U |: 1.9600 
Alfa: 5% 


U<|U|=> dados homogeneos 


**Teste Grubbs & Beck: 
*Limite superior: 

LS: 2458.6931 

Max: 1785.0000 

Alfa: 5% 


Max < LS => não há outlier 


*Limite inferior: 
LI: 183.2389 
Min: 246.0000 
Alfa: 5% 


Min > LI => não há outlier 
**Teste da assimetria: 
*Limite superior: 


Não há outlier - a retirada do máximo amostral não resultou em uma mudança significativa da assimetria. 


*Limite inferior: 


Não há outlier - a retirada do mínimo amostral não resultou em uma mudança significativa da assimetria. 


ESTIMAÇÃO DOS PARÂMETROS: 


Distribuição Posição Escala Forma 
Normal 746.5532  342.5656 
Log-Normal 6.5091 0.4801 

Gumbel 585.6067  278.8325 


Exponecial 360.0093  386.5439 
Pearson-IIl 746.5532  355.6089 1.0955 
Log-Pearson-lI | 6.5091 0.4810 -0.2474 
GEV 583.5028 2744976 -0.0167 
GPA 285.0576 | 640.4654 0.3878 
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TESTES DE AJUSTE: 

**Teste de Filliben 

Distribuição R(90%) R 
Normal 0.9803 0.9655 
Log-Normal 0.9803 0.9933 
Gumbel 0.9704 0.9946 
Exponecial 0.9641 0.9810 
Pearson-Il 0.9816 0.9949 
Log-Pearson-III 0.9850 0.9936 


GEV 
GPA 


0.9806 0.9948 
0.9867 0.9879 


**Teste da Variância 


Distribuição |Z| 
Normal 2.7931 
Log-Normal 0.6232 
Gumbel 0.1498 
Exponencial 2.0676 
Pearson-Il 0.3719 
Log-Pearson-llI 0.1078 
GEV 0.0149 
GPA 1.9071 
ESTIMAÇÃO DE QUANTIS: 


| Distribuição 
Normal 
Log-Normal 
Gumbel 
Exponencial 
Pearson-lllI 


Log-Pearson-lIl 


GPA 


2 


674.334 


5 10 50 100 


1051.827 1260.369 1574.320 1659.706 


200 


1724.967 


1805.161 
2958.870 
2511.574 
3030.160 
2406.210 
2509.553 
2593.246 
1823.213 
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BE CAARQUIV-AISEAFCLIPSCMD. EXE 
SISTEMA ESPECIALISTA 


1)Selecionando uma distribuição pelo teste da Variância dos momentos-L... 
a distribuição LogNormal é selecionada com 8.92 de confiança 
distribuição Gumbel é selecionada com 8.99 de confiança 
distribuição Pearson-III é selecionada com 6.97 de confiança 
distribuição LogPearson-III é selecionada com 6.99 de confiança 
distribuição GEU é selecionada com 1.88 de confiança 


2)Selecionando uma distribuição pelo teste de Filliben... 

a distribuição LogNormal é selecionada com 6.93 de confiança 

a distribuição Gumbel é selecionada com 6.98 de confiança 
distribuição Exponencial é selecionada com 8.86 de confiança 
distribuição Pearson-III é selecionada com 6.96 de confiança 
distribuição LogPearson-III é selecionada com 8.91 de confiança 
distribuição GEU é selecionada com 6.97 de confiança 
distribuição GPA é selecionada com 6.66 de confiança 


3)Procurando razões para rejeitar alguma das distribuições selecionadas... 
- a distribuição GPA é rejeitada, porque o parâmetro de forma 
é positivo 
- a distribuição Exponencial é rejeitada, porque o mínimo 
amostral é significativamente menor que o mínino teórico 
- a distribuição LogPearson-III é rejeitada, porque a assimetria dos 
logaritmos dos dados é negativa 


t4:01535 


E Fechar | Salvar | Dl [>] 
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AARQUIV- 1ISEAFACLIPSCMD. EXE 


a distribuição Pearson-III é selecionada com 9.96 de confiança 

a distribuição LogPearson-III é selecionada com 8.91 de confiança 
a distribuição GEU é selecionada com 6.97 de confiança 

a distribuição GPA é selecionada com 6.66 de confiança 


3)Procurando razões para rejeitar alguma das distribuições selecionadas... 
- a distribuição GPA é rejeitada, porque o parâmetro de forma 
é positivo 
- a distribuição Exponencial é rejeitada, porque o mínimo 
amostral é significativamente menor que o mínino teórico 
- a distribuição LogPearson-III é rejeitada, porque a assimetria dos 
logaritmos dos dados é negativa 


u)Verificando a parsimônia entre as distribuições selecionadas... 
- O aumento da incerteza na estimativa de mais um parâmetro 
para distribuição GEU não é justificado pelo teste de 
parsimônia em relação a distribuição Gumbel. 


S)Distribuições recomendadas... 
- Pearson-III com 8.9665 de confiança 
- Gumbel com 6.9837 de confiança 
- LogNormal com 6.9244 de confiança 


ATENÇÃO: feche esta janela antes de fazer outra análise. 
LIPS> À 


t4:01535 
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15 Iniciar É Docume: rdPad Sé Cia ISEAFL, e + 235 


Pelo relatório final apresentado, ver item 5, as três melhores distribuições são: Pearson II, Gumbel e Log- 
Normal, todavia a distribuição Gumbel apresentou coeficiente de ajuste mais próximo a 1, sendo a 
escolhida para fazer previsões futura. 


Como pode-se observar a conclusão foi praticamente a mesma, entretanto o resultado é obtido de forma 
rápida e prática. 
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5.1. OBJETIVOS INSTRUCIONAIS 
Ao final desta aula prática de Processamento Digital de Imagens você estará instruído para: 


1) Enunciar dois possíveis métodos aplicáveis à detecção de mudanças a partir de série temporal de 
imagens orbitais; 


2) Conceituar detecção de mudanças, NDVI (Índice de Vegetação por Diferença Normalizada) e RCEN 
(Radiometric Rotation Controlled by Nonchange axis) modificada e; 


3) Comparar os resultados obtidos pelos diferentes métodos de detecção de mudanças. 


5. 2.DETECÇÃO DE MUDANÇAS NO CONTEXTO DE SENSORIAMENTO REMOTO 


Realizar detecções de mudanças de características da superfície da Terra é importante para a 
compreensão tanto da dinâmica dos fenômenos quanto para a previsão dos impactos, bem como para o 
apoio na tomada de decisões. 


Durante as últimas décadas foram desenvolvidas várias técnicas de detecção de mudanças a partir de 
imagens, dentre elas as baseadas em imagens de Sensoriamento Remoto. Em geral, a detecção de mudança 
envolve a utilização de um conjunto de dados multitemporais que permite a análise qualitativa e 
quantitativa do fenômeno em estudo (CALDEIRA, 2016). 


Na prática, não é fácil selecionar um algoritmo adequado para um projeto de detecção de mudança, e por 
esta razão o conhecimento das técnicas de detecção de mudanças utilizados em pesquisas e aplicações 
anteriores é útil para entender como essas técnicas podem ser melhores exploradas e utilizadas na 
resolução de outros problemas. É relevante destacar que é necessário que estas mudanças sejam 
separáveis de outros fatores, como diferenças provocadas por condições atmosféricas, iluminação da cena, 
ângulos de visada do sensor, umidade do solo dentre outros. Visando isto, esta aula apresenta resultados 
da avaliação preliminar de duas abordagens de detecção de alterações. A primeira delas é feita a partir de 
uma modificação na técnica RCEN e a segunda baseada na diferença temporal entre imagens utilizando um 
índice de vegetação expressas pelo NDVI. Os dois métodos serão implementados e aplicados a um 
conjunto de imagens ortorretificadas, as quais foram tomadas em duas épocas distintas, ou seja, utilizando 
de série temporal. 


5.3.TÉCNICAS DE DETECÇÃO DE MUDANÇAS 


Como foi apresentado em sala de aula, há diversos métodos de detecção de mudanças, cada um com 
condições de aplicação específico. Vários autores resumiram algumas das técnicas mais utilizadas, com 
suas limitações e aplicações (SINGH, 1989; MOUAT et al. 1993; DEER, 1995; COPPIN e BAUER, 1996; 
JENSEN et al., 1997; SERPICO e BRUZZONE, 1999; LU et al., 2004, COPPIN et al., 2004; RADKE et al., 2005; 
THONFELD et al. 2005; WARNER et al. 2008; HUSSAIN et al., 2013; TEWKESBURY et al., 2015; MISHRA et 
al. 2017, ASOKAN e ANITHA, 2019). 


O objetivo da detecção de alterações é comparar a representação espacial de dois pontos no tempo, 
controlando todas as variações causadas por diferenças nas variáveis que não são de interesse, a fim de 
monitorar as mudanças causadas por diferenças nas variáveis de interesse (GREEN; KEMPKA e LACKEY, 
1994). 


A premissa básica da utilização de dados de Sensoriamento Remoto na detecção de alterações é que as 
mudanças nos objetos de interesse resultam em variações nos valores de reflectância ou na textura local, 
mas essas análises podem ser separadas de mudanças causadas por outros fatores como diferenças nas 
condições atmosféricas, iluminação, ângulo de visada e umidade do solo (DEER, 1995). 


Em sala de aula foi visto que alguns processamentos prévios são necessários para que a detecção de 
mudanças seja realizada de forma correta, são eles: Correção Geométrica, Radiométrica e Atmosférica. 


Como o foco da nossa aula prática é a aplicação e comparação entre duas técnicas, já aplicado às correções 
necessárias previamente as imagens, serão considerados as técnicas descritas a seguir. 
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5.3.1. MÉTODO DE DETECÇÃO DE MUDANÇAS BASEADO NA DIFERENÇA TEMPORAL DE NDVI 


Z 


Para viabilizar a detecção de mudanças, segundo essa abordagem, é necessário realizar uma prévia 
normalização radiométrica nas bandas espectrais e imagens multitemporais envolvidas no processo, além 
da correção geométrica entre as bandas. 


O Índice de Vegetação por Diferença Normalizada (NDVI) conforme (LEE et al. 2020), pode ser obtido a 
partir da normalização do índice de razão simples a partir da Equação (1). Ressalta-se, ainda, que o 
resultado dessa aplicação do NDVI será representado entre um intervalo de -1 a 1, ou seja, as análises que 
serão relacionadas aos parâmetros do NDVI deverão estar enquadradas neste intervalo. Os fatores que 
corresponderão a esses dois índices numéricos de representação estarão diretamente ligados ao 
quantitativo de vegetação presente na cena, uma vez que, segundo aborda LIU (2007), quanto mais densa 
estiver essa vegetação o valor se amostrará mais próximo de 1. Quanto ao contrário, ou seja, índices mais 
próximos de 0, ter-se-ia esse valor quando não se apresenta uma presença expressiva de vegetação sendo, 
assim, representados por solo exposto, construções, áreas descampadas dentre outros. Tratando-se de 
valores negativos, ou seja, mais próximos de -1 estaria sendo representando corpos d'água e sombra. 


Pnir7 Pred 
NDVI = rtdar trtiteos (1) 
PnirtPred 
onde, pnir é o fluxo radiante refletido no infravermelho próximo e prea é o fluxo radiante refletido no 


vermelho da região do visível. 


Com base nos valores de NDVI calculados para cada sensor é calculado o valor que deve indicar possíveis 
alterações, como visto na Equação (2). 


NDVIetecção = NDVlIp, — NDVIp, (2) 


sendo: NDVlIetecção à imagem NDVI de detecção de alterações; NDVIp; e NDVIpz os valores de ND das 
imagens NDVI de cada uma das datas. 


5.3.2. TÉCNICA RCEN (RADIOMETRIC ROTATION CONTROLLED BY NONCHANGED AXIS) 
MODIFICADA 


Essa categoria tem como referência a técnica de detecção de mudanças RCEN, a qual faz uma 
transformação representando uma similaridade entre a resposta espectral dos pixels de não mudança nas 
épocas em que se adotaram as coletas. A técnica de origem foi proposta por (MALDONADO; SANTOS; 
GRAÇA, 2007), e desse modo, aplicando-se no método original, as seguintes premissas são verdadeiras: 


= A maior parte dos pixels entre dois rasters, presente em uma mesma região na banda espectral, 
em um curto período de tempo, não é alterado, principalmente se tratando de regiões monitoradas; 


= O eixo que representa a não mudança gera um ângulo 6 indicando uma posição relativa que não 
sofre alteração; 


= O eixo que representa a não mudança pode ser obtido por uma regressão linear de valores 
radiométricos, os quais são coletados em locais onde não sofreram mudanças; 


n As bandas espectrais corretas devem ser selecionadas conforme o alvo de interesse; 
= A compatibilidade radiométrica entre as bandas envolvidas é garantida; 
= A técnica é robusta até quando se compara dados de sensores diferentes. 


Em concordância com os pressupostos destacados por (MALDONADO; SANTOS; GRAÇA, 2007), o 
resultado do processo de aplicação da técnica será o eixo de não mudança, sendo este obtido por meio de 
informações nas áreas onde não sofreram alterações. Desse modo, sendo assumido que o nível de 
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alterações não é o predominante, ao se realizar a comparação das imagens dessas duas épocas, é possível 
conceituar que o eixo de não mudança pode ser de forma estimada a partir do raster de forma geral, ou 
seja, sem que seja realizada uma pré-seleção das áreas sem modificações. 


Dessa forma, descreve-se o método de detecção denominado como RCEN modificado. O ângulo de rotação 
8 (calculado pelo eixo que representa a não mudança) pode ser obtido se aplicando uma regressão linear 
nos valores de ND (Nível Digital) das bandas espectrais correspondentes, sendo que variações mínimas na 
largura da banda ou na sua sensibilidade podem ser assimiladas pela transformação aplicada. Deste modo, 
tendo-se como sabe duas imagens em diferentes épocas, tem-se a Figura (1), a qual traz representações de 
os elementos essenciais para o desenvolvimento da equação que foi calculada. 


Figura 1. Representação gráfica do eixo de não mudança e do valor de Ir (imagem intermediária) 
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Fonte: (CALDEIRA et al. 2018) 


Ao se considerar um pixel com valores determinados por Ip, e Ip, nos dois períodos - respectivamente, 
se este pixel estiver situado sobre a reta representada pelo coeficiente angular tan(6), isso nos direciona a 
interpretar que esse pixel não sofreu alterações. Se esse pixel estiver disposto como representado na 
Figura (1) (ponto em vermelho), a distância entre as projeções (representado em azul na figura) será um 
indicador de que este pixel sofreu modificações entre as épocas escolhidas. 


Ou seja, o valor de Ir, presente na Equação (3), será estimado por uma diferença entre as distâncias hz - 
hy. Dessa forma, ao se aplicar as funções trigonométricas adequadas para cada um dos triângulos, é viável 
se determinar os valores de h; e h; em função dos valores de In1, Ipz, 8, bem como de b (coeficiente linear 
da reta de não mudança). Ou seja: 


Inr = (cos0)"!+ Ip — (sin60)!+ (Ib2—b) 3) 


A utilização do coeficiente linear (b), presente na Equação (3), pode ser justificada pela razão de que, na 
Figura (1), a reta passa pela origem, ou seja, apresentando um caso ideal. Se a reta não passar pela origem, 
será necessário realizar uma translação desta reta para a origem utilizando o coeficiente linear. Essa 
operação resulta em uma imagem intermediária (Iwnr) a qual se relaciona com as modificações ocorridas. 
Logo, quanto maior for a mudança, maior será o valor resultante. 


Desse modo, pode-se notar, na Figura (1) conforme apresentada, que os pixels analisados podem ficar 
acima ou abaixo desta reta formada pelo coeficiente angular tan(0), resultando em valores os quais podem 
ser positivos e negativos para Ir. 
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Baseando-se em valores intermediários que foram calculados para cada banda espectral utilizadas nas 
duas épocas, em base da Equação (3), calcula-se o valor indicando possíveis alterações. Na Equação (4), 
pode-se obter o resultado com o índice de alteração Ietec calculado a partir de das bandas multiespectrais: 


ps ist 
Lo Da Irei (4) 


sendo Irei Os valores intermediários para as bandas; i os índices das bandas e n o número de bandas. 


Os valores obtidos pela Equação (3), para uma única banda, resultarão em um número pertencente ao 
conjunto dos reais. Consonante a isso, o resultado desses modelos matemáticos pode ser de variação 
positiva à negativa, o qual poderá indicar o aparecimento de “novos” alvos bem como o “desaparecimento” 
de alguns objetos na cobertura da terra. Para realçar as mudanças mais expressivas, pode ser adotado de 
um percentual, ou seja, de um limiar, o qual permitirá esse realce às mudanças de interesse. 


5.4. PRÁTICA 


Nossa aula prática será realizada em Laboratório, em virtude de sua estrutura a qual possui computadores 
e softwares específicos para a atividade. 


5.4.1. MATERIAIS 
Em síntese, os materiais que serão utilizados para esta aula prática serão: 
= Imagens ortorretificadas orbitais do mesmo local, porém em diferentes épocas; e 


= Um computador desktop. 


5.4.2. ANTES DA MEDIÇÃO 


Primeiramente, em aulas passadas, os pares de cenas do mesmo sensor, referentes às bandas 
correspondentes nas duas datas, foram co-registradas geometricamente utilizando o software ENVI 5.0 e 
realizado a correção atmosférica antes de serem submetidas ao processamento de análise de alterações. 
Lembrando que a detecção de mudanças, utilizando as técnicas descritas, faz a comparação pixel a pixel, 
logo o tamanho do pixel terá que ser o mesmo para as duas imagens. Caso as imagens utilizadas forem de 
sensores diferentes, além do registro geométrico entre as imagens, também deverá haver uma 
reamostragem na imagem de maior resolução espacial de modo a compatibilizá-las. Um exemplo são os 
sensores SPOT-6, com GSD (Ground Sample Distance) igual a 6 metros e o WorldView-2 (GSD de 2 metros). 
A Figura (2) mostra um fluxograma do método proposto para se promover uma melhor compreensão dos 
procedimentos. 
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Figura 2. Fluxograma do método proposto 
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erecorte-das 
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As imagens possuem o g Reamostrar a imagem 
mesmo tamanho:de pixelZ de maior resolução 


Aplicar um:método de 
detecção de mudança 


Aplicação de 
um limiar para 
orealcedas 
mudanças 


Exibir as 
detecções 


Fonte: Autores (2020) 


Os segmentos de imagem devem possuir o mesmo tamanho, com a mesma área geográfica e mesmo 
tamanho do pixel, porém em diferentes datas para cada sistema imageador. 


As Figuras (3) e (4) apresentam exemplos de recortes de imagens que representam áreas testes, na forma 
de composições coloridas RGB. As imagens indicadas pelas letras (a) e (b) em cada figura, diferenciam as 
cenas dos sensores WorldView-2 e Spot-6, respectivamente, nas quais se percebe com clareza o melhor 
detalhamento geométrico da imagem WorldView-2, definido pelo menor GSD desse sensor. 


Figura 3. Área teste extraída de composição colorida normal (RGB) das imagens: (a) WordView-2 de 
29/07/2010; (b) Spot-6 de 22/04/2013. 


(a) (b) 
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Fonte: Autores (2020) 


Figura 4. Área teste extraída de composição colorida normal (RGB) das imagens: (a) WordView-2 de 
02/09/2011; (b) Spot-6 de 22/10/2012. 


(a) (b) 


Nota-se que nas Figuras (3) e (4) há a presença de alterações entre as duas datas, como: trapiches, casas e 
áreas de exploração irregular do solo, as quais devem ser detectadas pelo algoritmo de detecção de 
alterações. 


Para a visualização das mudanças mais significativas, tanto as que correspondem a novos objetos /feições 
quanto a objetos/feições que deixaram de existir, serão aplicados limiares para realçar as alterações 
detectadas. Estes limiares são valores aleatórios resultantes das imagens de detecção quando o usuário 
sabe onde há uma alteração. Nestas situações, ele pode usar este valor do pixel como base. Inicialmente 
define-se um valor percentual Lperc que corresponde às maiores modificações. Assumindo-se que as 
imagens que permitem detectar alterações, seja pela Equação (2) ou Equação (4), estejam normalizadas 
no intervalo [0-255], a partir do valor Lperc pode-se determinar os seguintes limiares superiores e 
inferiores (expressos em valores de ND), com a Equação (5): 


100 — a 
100 


Line = 255 — Loyp 


Lsup — 255 * ( 
(5) 


sendo: Lsype Linr OS limiares superior e inferior. Para efeito de visualização, aos valores de ND maiores que 
Lsupserão realçados na cor azul e os valores de ND menores que L;yr foram associados à cor vermelha. 


Para o cálculo das imagens e aplicação das equações, pode-se utilizar qualquer linguagem de programação, 
como Pyhton, Matlab, C, C++, Java, entre outras. A metodologia é simples, o usuário terá que fazer com que 
o programa leia as imagens já corrigidas previamente e após isso executar, pixel a pixel, as equações 
mostradas anteriormente. 


A Figura (5) exemplifica um código em Matlab para a aplicação da Equação (1), onde as variáveis i e j 
representam as posições dos pixels em linha e coluna, respectivamente e le c, o total de linhas e colunas 
das imagens, respectivamente. 
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Figura 5. Exemplo de parte de um código em Matlab, para o cálculo do NDVI. 


IV = double (imagem(:,:,4));% Lendo a banda do Infravermelho Proximo 
= double (imagem(:,:,3));% Lendo a banda do Vermelho 


[1,c] = size(IV);%Retirando o tamanho da imagem 


for i=l:1 
for j=":c 
(1,5) = (IVv(i,5)-V(i,5))/(Iv(i,5)+V(i,5));S%Aplicando a Equação do NDVI para cada pixel 


Fonte: Autores (2020) 


As Figuras (6) e (7) apresentam as imagens NDVI (Equação 1) geradas para cada área teste e para ambos 
os sensores utilizando o código da Figura (5). 


Figura 6. Imagens NDVI da área teste resultantes da aplicação da Equação (1) às imagens multissensor: (a) 
Wordview-2, (b) Spot-6. 


(a) (b) 
Fonte: Autores (2020) 


Figura 7. Imagens NDVI da área teste resultantes da aplicação da Equação (1) às imagens multissensor: (a) 
Wordview-2, (b) Spot-6. 


A 


Fonte: Autores (2020). 
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5.4.3. COMPARAÇÃO DOS RESULTADOS A PARTIR DA APLICAÇÃO DE LIMIARES 


Para ter uma melhor interpretação dos dados e consequentemente uma melhor avaliação do algoritmo, o 
usuário deverá utilizar de um limiar para destacar as alterações que foram construídas de um ano para 
outro, ou seja, alterações com grandes diferenças positivas nos ND's e alterações que foram removidas de 
um ano para o outro, ou seja, alterações com grandes diferenças negativas nos ND's. 


As Figuras (8) e (9) apresentam as detecções realçadas em ambos os métodos para as áreas testes, os 
quais foram sequenciados a partir de uma imagem representando a detecção principal a ser encontrada 
pelos métodos citados. Estes limiares não são válidos para representar a validação dos métodos. O 
procedimento de validação será realizado em uma próxima etapa, assim sendo, estes limiares serão 
apresentados somente para uma melhor visualização dos resultados, como dito anteriormente. 


Figura 8. Imagem da área teste com detecções realçadas, sendo em vermelho as alterações que foram 
removidas ao longo do tempo e em azul as alterações construídas ao longo do tempo, onde em (a) imagem 
final do método RCEN modificado, com destaque dos alvos a serem detectados, em (b) imagem final do 
método de diferenciação de imagens NDVI, com destaque dos alvos a serem detectados e em (c.1 e c.2) 
imagens exemplo, para os alvos a serem detectados. 


(c.2) 


Fonte: Autores (2020) 


A Figura (8) mostra uma área onde se destacam dois grupos de alterações. As alterações de áreas com 
construções e áreas com remoções, utilizando-se de um limiar de 30%, ou seja, os 30% maiores valores de 
ND foram destacados. Nesta área destacam-se as construções representadas na Figura (8.c.1) que foram 
removidas entre os anos de 2011 e 2014 e na Figura (8.c.2) que foram construídas. Respeitando a lógica 
anterior, as alterações teriam que aparecer em tons de vermelho e azul, respectivamente. Percebe-se que 
utilizando o método NDVI os alvos construídos não foram destacados pontualmente e nem corretamente, 
porém este método destacou alguns alvos que foram removidos. O fato de ter crescido vegetação no lugar 
dos alvos removidos, ajudou para detectar estes alvos. Comparando com o método RCEN modificado 
percebe-se que este destacou os alvos construídos corretamente, ou seja, da cor azul e bem definido na 
imagem, representando então que este alvo não estava presente em 2011 e agora está em 2014. 
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Figura 9. Imagem da área teste com detecções realçadas, sendo em vermelho as alterações que foram 
removidas ao longo do tempo e em azul as alterações construídas ao longo do tempo, onde em (a) imagem 
final do método RCEN modificado, com destaque dos alvos a serem detectados, em (b) imagem final do 
método de diferenciação de imagens NDVI, com destaque dos alvos a serem detectados e em (c) imagem 
exemplo para os alvos a serem detectados. 


(a) (b) 
Fonte: Autores (2020) 


A Figura (9) mostra uma área onde se destacam duas áreas de alterações, utilizando-se de um limiar de 
30%, ou seja, serão destacados os 30% menores valores de ND, onde as alterações representam uma 
extração ilegal de argila. Percebe-se que ambos os métodos obtiveram bons resultados, pois detectaram 
boa parte da extração de argila. Porém nota-se que o método RCEN modificado ainda demonstra melhores 
resultados em relação ao método de NDVI. 


Este método de destaque das detecções pode ser feito apenas modificando o valor original do pixel pelas 
cores vermelha ou azul, ou seja, basta apenas percorrer a imagem e quando o valor do pixel for maior ou 
igual ao limiar proposto, adota-se as cores azul ou vermelha, como visto em um exemplo de código na 
Figura (10). 


Após a execução dos códigos, discute-se as vantagens e desvantagens de cada método, como no exemplo 
mostrado. 


Figura 10. Exemplo de parte de um código em Matlab para destaque em vermelho e azul para os pixeis que 
apresentaram alteração 


[lin,col] = size (imagemDIF);SVerificando o tamanho da imagem 
Limiar = 255.*(100. - Lperc)/100;%calculo do limiar superior 
Limiar2 = 100.*(100. - Lperc)/100;%calculo do miliar inferior 
imagemBIN=zeros (lin,col,3); 
:,:,2) = imagemBASE(:,:); 
:,:,3) = imagemBASE(:,:); 


for j=!:col 
if imagemDIF(i,5j) >= LimiarSAdotando a cor Vermelha 
imagemBIN(i,5,1)=0; 
imagemBIN(i,5,2)=0; 
imagemBIN(i,5,3)=255; 
else 
if imagemDIF(i,5j) <= Limiar2%Adotando a cor azul 
imagemBIN(i,5j,1)=255; 
imagemBIN(i,5,2)=0; 
imagemBIN(i,5,3)=0; 
end 
end 
end 
end 


Fonte: Autores (2020). 


Práticas em Saneamento, Clima e Meio Ambiente 


BIBLIOGRAFIA 


[1] ASOKAN, A.; ANITHA, J. Change detection techniques for remote sensing applications: a survey. Earth Science 
Informatics, /S. LJ, n. 12, p. 143-160, 2019. DOI: 10.1007/s12145-019-00380-5. 


[2] CALDEIRA, C. R. T. Detecção de mudanças no entorno de reservatórios a partir de série temporal de imagens 
orbitais. 2016. Dissertação (Mestradoem Ciências Cartográficas) - - FCT/Universidade Estadual Paulista, Presidente 
Prudente, 2016. 


[3] CALDEIRA, C. R. T.; GALO, M.; IMAI, N. N.; GALO, M. de L. T. B.; HASEGAWA,)J. K.; AMORIM, A.; SHIMABUKURO, 
M. H.; PINTO, M. Detecção de mudanças de forma automática no entorno de reservatórios: estudo de caso - Canoas 1. 
Revista Brasileira de Cartografia, [S. LJ, v. 70, n. 4, p. 1348-1377, 2018. DOI: 10.14393/rbcv70n4-46508. 


[4] COPPIN, P.; JONCKHEERE, I.; NACKAERTS, K.; MUYS, B.; LAMBIN, E. Digital change detection methods in 
ecosystem monitoring: A review. International Journal of Remote Sensing, /S. LJ, v. 25, n. 9, p. 1565-1596, 2004. DOI: 
10.1080/0143116031000101675. 


[5] COPPIN, P. R.; BAUER, M. E. Digital Change Detection in Forest Ecosystems with Remote Sensing Imagery. 
Remote Sensing Reviews, /S. L], v. 13,n. 3-4, p. 207-234, 1996. DOI: 10.1080/02757259609532305. 


[6] DEER, P. J. Digital Change Detection Techniques: Civilan and Military Applications. In: C 
Org. INTERNATIONAL SYMPOSIUM ON SPECTRAL SENSING RESEARCH 1995, Report (Greenbelt, MD: Goddard Space 
Flight Center). Anais [...]. Report (Greenbelt, MD: Goddard Space Flight Center) p. 1-13. 


[7] GREEN, K.; KEMPKA, D.; LACKEY, L. Using remote sensing to detect and monitor land-cover and land-use 
change. Photogrammetric Engineering and Remote Sensing, /S. L.], v. 60, n. 3, p. 331-337, 1994. 


[8] HUSSAIN, M.; CHEN, D.; CHENG, A.; WEI, H.; STANLEY, D. Change detection from remotly sensed imagens: 
From pixel-based to object-based approaches. ISPRS Journal of Photogrammetry and Remote Sensing, /S. LJ, v. 80, p. 
91-106, 2013. 


[9] JENSEN, J. R.; COWEN, D.; NARUMALANI, S.; HALLS, J. Principles of change detection using digital remote 
sensor data. In: (. Org. JINTEGRATION OF GEOGRAPHIC INFORMATION SYSTEMS AND REMOTE SENSING 1997, 
Cambridge. Anais [...). Cambridge: Cambridge University Press, 1997. p. 37-54. 


[10] LEE, K.; KIM, K.; LEE, S. G.; KIM, Y. Determination of the normalized difference vegetation index (NDVI) with 
top-of-canopy (TOC) reflectance from a KOMPSAT-3A image using orfeo tool box (OTB) extension. ISPRS International 
Journal of Geo-Information, /S. L], v. 9,n. 4, p. 1-16, 2020. DOI: 10.3390/ijgi9040257. 


[14] LIU, W. T. H. Aplicações de sensoriamento remoto. UNIDERP ed. Campo Grande: UNIDERP, 2007. v. 53 DOI: 
10.1017/CB09781107415324.004. 


[12] LU, D.; MAUSEL, P.; BRONDIZIO, E.; MORAN, E. F. Change Detection Techniques. International Journal of 
Remote Sensing, /S. .], v. 25, n. 12, p. 2365-2407, 2004. 


[13] MALDONADO, F. D.; SANTOS, J. R.; GRAÇA, P. M. L. Change detection technique based on the radiometric 
rotation controlled by no-change axis, applied on a semi-arid landscape. International Journal of Remote Sensing, [S. 
1], v. 28,n.8, p. 1789-1804, 2007. DOI: 10.1080/01431160600905821. 


[14] MISHRA, S.; SHRIVASTAVA, P.; DHURVEY, P. Change Detection Techniques in Remote Sensing: A Review. 
International Journal of Wireless and Mobile Communication for Industrial Systems, [S. LJ, v. 4,n. 1, p. 1-8, 2017. DOI: 
10.21742 /ijwmcis.2017.4.1.01. 


[15] MOUAT, D. A.; MAHIN, G. G.; LANCASTER, J. Remote sensing techniques in the analysis of change detection. 
Geocarto International, [S. LJ, v. 8, n. 2, p. 39-50, 1993. DOI: 10.1080/10106049309354407. 


[16] RADKE, R.).; ANDRA, S.; AL-KOFAHI, O.; ROYSAM, B. Image change detection algorithms: A systematic survey. 
IEEE Transactions on Image Processing, [S. LJ, v. 14, n. 3, p. 294-307, 2005. DOI: 10.1109/TIP.2004.838698. 


[17] SERPICO, S. B.; BRUZZONE, L. Change detection. In: (., Org. JINFORMATION PROCESSING FOR REMOTE 
SENSING 1999, Singapore: World Scientific Publishing. Anais [...]. Singapore: World Scientific Publishing p. 319-336. 


[18] SINGH, A. Digital change detection techniques using remotely-sensed data. International Journal of Remote 
Sensing, /S. 1.], v. 10, n. 6, p. 989-1003, 1989. 


[19] TEWKESBURY, A. P.; COMBER, A.J.; TATE, N.).; LAMB, A.; FISHER, P. F. A critical synthesis of remotely sensed 
optical image change detection techniques. Remote Sensing of Environment, [S. LJ], v. 160, p. 1-14, 2015. DOI: 
10.1016/j.rse.2015.01.006. 


[20] THONFELD, F.; HECHELTJEN, A.; BRAUN, M.; MENZ, G. From algorithms to processing chains: A review of 
land cover and land use change detection methodologies. Technology, /[S. L.], n. 1, p. 99775-99775, 2005. 


[21] WARNER, T. A.; ALMUTAIRI, A.; LEE, J. Y. Remote Sensing of Land Cover Change. In: (. Org.)THE SAGE 
HANDBOOK OF REMOTE SENSING 2008, Anais [...]. [s.l: s.n.] p. 459-474. DOI: 10.4135/9780857021052.n33. 


Capítulo 6 


Aula prática: Nivelamento geométrico, 
determinação de desníveis e altitudes 


Mayara Cobacho Ortega Caldeira 
Carlos Rodrigo Tanajura Caldeira 
Lucas Daniel Noronha Ferreira 


Práticas em Saneamento, Clima e Meio Ambiente 


6.1. OBJETIVOS INSTRUCIONAIS 


Ao final deste capítulo teórico e prático referente ao Nivelamento Geométrico, determinação de desníveis 
e de altitudes, os alunos serão capazes de: 


4) Enunciar este método aplicável à observação de altitudes do terreno; 

5) Identificar o uso de equipamentos necessários para a execução; 

6) Conceituar Nivelamento Geométrico, desnível e altitude, e; 

7) Operacionalizar as leituras em campo e cálculos a partir do emprego deste método. 


6.2. NIVELAMENTO GEOMÉTRICO 


A determinação de um ponto da superfície terrestre só está completa, no ponto de vista da Geomática, se 
além das coordenadas planimétricas (X, Y), estiver também associada a elas uma componente altimétrica 
(h ou H) (CORDINI, 2014). Conforme afirma Silva e Segantine (2015), é por meio da altimetria que se 
representa o relevo de um terreno para projetos de engenharia e mapeamento, como por exemplo 
projetos de vias de transporte, saneamento, distribuição de água e esgoto, construção de barragens, entre 
outros. 


A determinação de desníveis entre dois pontos e o transporte de altitudes em qualquer parte da superfície 
terrestre representa uma das operações fundamentais sob atuação da Topografia e Geodésia, conforme 
afirma Ghilani e Wolf (2012). Atualmente, com o surgimento de novas tecnologias, houve um aumento na 
diversidade de se obter esse desnível acompanhada de um aumento na quantidade e variedade de 
equipamentos, os quais cada um possui sua precisão que possibilitaram o desenvolvimento e 
aprimoramento da técnica (NICACIO et al. 2018). 


De acordo com Norma Brasileira NBR13.133, o conjunto de operações e técnicas realizadas para 
determinar a altitude, cotas, ou desníveis dos pontos de interesse, é denominada de Levantamento 
Altimétrico ou Nivelamento (ABNT, 1994). 


Conforme visto em outras aulas e destacado Silva e Segantine (2015) e em outras literaturas, existem 
vários métodos de nivelamento disponíveis, como por exemplo: Nivelamento Trigonométrico, 
Nivelamento Taqueométrico, Nivelamento Geométrico, Nivelamento por Satélite GNSS e Nivelamento a 
Laser e Digital, e também o método antigo de nivelamento Barométrico o qual se baseava na medição de 
pressão atmosférica associada à altitude. No entanto, o enfoque deste capítulo é o Nivelamento 
Geométrico. 


O Nivelamento Geométrico é definido como o método de determinação das componentes altimétricas com 
base nos dados obtidos através de visadas horizontais a partir da diferença de nível entre pontos do 
terreno (GHILANI; WOLF, 2012). Conforme apresentado Figura 1, este processo é realizado por 
intermédio de leituras em miras graduadas colocadas verticalmente nos referidos pontos de interesse, 
obtidas com equipamento topográfico denominado Nível (ABNT, 1994; VEIGA; ZANETTI; FAGGION, 2012). 


Figura 1. Princípio Básico do nivelamento Geométrico. 


Plano de Referência 


Horizontal 


Fonte: Adaptado de Silva e Segantine (2015) 
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O Nível deve ser instalado e nivelado para que seu eixo vertical coincida com a vertical do lugar, deste 
modo, a luneta estabelece um plano horizontal de referência (KAVANAGH; MASTIN, 2014). Vale ressaltar 
que diferentemente dos outros métodos de nivelamento, o desnível é determinado com base em ângulos 
verticais ou zenitais. Além disso, neste método as distâncias verticais (desnível) são obtidas a partir da 
interseção do plano horizontal com a mira vertical (SILVA; SEGANTINE, 2015). 


Deste modo, o nivelamento geométrico estabelece operações de obtenções dos desníveis em pontos de 
interesse no terreno com suas alturas referidas a uma superfície arbitrária (cotas) ou ao nível médio do 
mar (altitude) (ABNT, 1994). 


Diante disso, a relação que é atribuída ao uso deste equipamento durante o levantamento possibilita a 
determinação do desnível entre dois pontos, podendo ser aplicada em duas vertentes diferentes, sendo 
elas: 1) geodésicas ou 2) topográficas (VEIGA; ZANETTI; FAGGION, 2012). Porém, convém ressaltar que 
dentro da área de levantamentos geodésicos este método de levantamento é utilizado para fins 
monitoramentos de grandes estruturas, o que embasa o necessário entendimento acerca desta forma de 
nivelamento (SIGUEL et al., 2013). 


6.3. EQUIPAMENTOS: NÍVEL TOPOGRÁFICO E MIRA ESTADIMÉTRICA 


Conforme visto da Figura 2, o Nível é um equipamento que possui uma luneta telescópica e pode girar em 
torno do seu eixo vertical. Além disso, possui uma base nivelante com um ou dois níveis de bolha para ser 
fixado a um tripé (NADOLINETS, LEVIN; AKHMEDOV, 2017; SILVA; SEGANTINE, 2015). Neiga, Zanetti e 
Faggion (2012), definem o Nível como: 


“Equipamentos que permite definir com precisão um plano horizontal 
ortogonal à vertical definida pelo eixo principal do equipamento. As principais 
partes de um nível são: 1) Luneta; 2) Nível de bolha; 3) Sistemas de 
compensação (para equipamentos automáticos); 4) Dispositivos de calagem” 
(VEIGA; ZANETTI; FAGGION, 2012, p. 196). 


Figura 2. Representação do Nível Topográfico 
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Fonte: Autores (2020) 


No entanto, atualmente - além destes níveis ópticos, com o advento da tecnologia pode-se empregar 
equipamentos mais robustos e precisos, como Níveis Digitais e Níveis Laser (NADOLINETS; LEVIN; 
AKHMEDOV, 2017). Todavia, a NBR 13.133 (ABNT, 1994) classifica os níveis quanto a precisão, 
considerando 1 km de duplo nivelamento, conforme apresentado na Tabela 1. 


Tabela 1. Classificação dos níveis segundo a NBR13.133 


Classes de Níveis Desvio-Padrão 


1- Precisão Baixa >+ 10 mm/km 
2- Precisão Média <+10mm/km 
3- Precisão Alta <+03 mm/km 
4- Precisão Muito Alta <+01mm/km 


Fonte: Adaptado (ABNT, 1994, p.6) 
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Outro equipamento que integra o nivelamento geométrico é a Mira estadimétrica, ou também denominada 
de Falante ou Régua estadimétrica. Porém, existem miras para operações normais (Figura 3b) e miras de 
precisão mais elevada (de ínvar, empregadas na rede altimétrica do Instituto Brasileiro de Geografia e 
Estatística- IBGE), ou ainda, miras com códigos de barra quando se utiliza níveis Digitais (Figura 3a) 
(CORDINI, 2014; KAVANAGH; MASTIN, 2014). 


Figura 3. Mira estadimétrica 


a A HOUTEN) CUCZUUO 


b [iii rossi retas 


Fonte: Furtado, Schmidt (2020) 


As miras têm a grande vantagem de possibilitar que sejam determinadas as leituras diretamente do 
aparelho com auxílio dos fios estadimétricos (Superior - FS, Médio - FM e Inferior - FI) do retículo da 
luneta, as alturas de visadas aos pontos topográficos (GARCIA; PIEDADE, 1984; VEIGA, ZANETTI; 
FAGGION, 2012), como apresentado na Figura 4a. 


Figura 4. Fios estadimétricos (a) e graduações da Mira (b 


Fonte: Autores (2020) 


As miras devem ser colocadas verticalmente nos pontos a nivelar, podendo ser de encaixe ou dobráveis. 
Vale ressaltar que para manter a verticalidade, um nível de cantoneira ou um nível de bolha junto à 
mesma, faz-se necessário (KAVANAGH; MASTIN, 2014). Sua menor célula gráfica de graduação é o 
centímetro e são numeradas de decímetro em decímetro (Figura 4b) (CORDINI, 2014). 


A leitura na mira é constituída de um número de quatro casas decimais (metro, decímetro, centímetro e 
milímetro estimado pelo operador), como o exemplo a Leitura do Fio Médio da Figura 4a é 1,431m, sendo 
que o milímetro é interpretativo e varia para cada operador (CORDINI, 2014). Para evitar erro, deve-se 
observar o sentido do crescimento da graduação. Além disso, frisa-se que um erro de 1 mm na leitura da 
mira ocasiona um erro de 10 cm no terreno (IBGE, 2017). 


6.4. CÁLCULOS MATEMÁTICOS PARA DETERMINAÇÃO DA ALTITUDE 


O procedimento de campo é realizado a cada lance, que é o conjunto básico formado pelas visadas a rée a 
vante, ou seja, as etapas envolvendo mira a mira (IBGE, 2017). Tais processos consistem em nivelar e 
instalar o Nível entre os pontos de interesse, visar a mira posicionada verticalmente no ponto inicial, o 
qual possui altitude conhecida que será a referência altimétrica para os demais pontos, e realizar a leitura 
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dos fios estadimétricos, denominada de leitura da ré ou visada à ré (SILVA; SEGANTINE, 2015). Em 
seguida, repete-se o procedimento de leitura para a vante. 


Considerando a Figura 5, é fácil notar que a diferença entre a Leitura da ré e a vante resulta no desnível 
(4h), conforme expresso na Equação 1. 


Aha — FMaRé = FMp't (1) 


Sendo: 
FM»: fio médio da visada à ré; 


FMg*t : fio médio da visada à vante. 


Figura 5. Lance do Nivelamento Geométrico 


Fonte: Adaptado de Silva e Segantine (2015) 


Conhecendo a altitude do ponto A (Ha) e considerando o desnível obtido na Equação 1, pode-se 
determinar o valor da altitude do ponto B (Hs) a partir do transporte de altitudes conforme o cálculo das 
altitudes prévias apresentadas na Equação (2). 


Hg = Ha + Aha (2) 


No entanto, há dois tipos de nivelamento geométrico: simples e composto. No nivelamento geométrico 
simples (Figura 6), a medição dos desníveis e altitudes dos pontos são determinadas a partir de uma única 
instalação do Nível, tomando como referência a leitura à ré para os demais pontos (GHILANI; WOLF, 2012; 
VEIGA, ZANETTI; FAGGION, 2012). Um limitante para este método é a distância entre os pontos (não 
podendo ser muito grande) e a variação do terreno, uma vez que os pontos devem ser visados de uma 
única estação, ou seja, aplicações se restringem a pequenas áreas e relativamente planas (SILVA; 
SEGANTINE, 2015). 
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Figura 6 - Nivelamento Geométrico Simples 


Fonte: Autores (2020) 


Já no que concerne ao nivelamento geométrico composto (Figura 7), estaciona-se o aparelho em mais de 
uma estação para se nivelar o local de estudo (CORDINI, 2014). Ou seja, é quando ocorre uma sucessão de 
nivelamentos por um caminhamento, os quais são realizadas leituras a ré de pontos intermediários, 
denominados de pontos de mudança (SILVA; SEGANTINE, 2015). 


Figura 7. Nivelamento Geométrico Composto 


A 


Sentido do Camnhamento 


Fonte: Autores (2020) 


Conforme indicado na Figura 7, tem-se a seguinte sequência de cálculos para determinar o desnível entre 
o Ae B (Equação 3). Nota-se que os pontos 1 e 2 são pontos de mudança, que serão referência a ré para 
determinação dos pontos 2 e B, respectivamente. 


Ahas =— FMAa- FM, — Hi = Ha + Ahas 
Ahy2 = FM,- FM, 5 H2=H4 + Ahy2 
Ah» — FMs- FMp > > Hp = H, + Ah» 


Aha = >» FMPé E » piputisto ou ainda, 
Aha,p = 5h (3) 
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Além disso, vale frisar que esta metodologia pode ser empregada para vários pontos de mudança, de modo 
a abranger toda área de estudo. Assim como, possibilita que sejam efetuadas irradiações em cada ponto 
estacionado, similar ao procedimento do nivelamento simples, afim de um detalhamento altimétrico da 
superfície a ser representada, conforme explicitado na Figura 8. Na literatura, este método é denominado 
de Caminhamento Misto ou Nivelamento com irradiações. 


Figura 8. Nivelamento Geométrico com Caminhamento Misto 


Fonte: Autores (2020) 


Na Figura 8, os pontos A, P1, P2, P3, P4, B pertencem ao caminhamento sendo que o ponto A e B devem ter 
altitude conhecida para possibilitar o erro de fechamento e transporte de altitude. Já os pontos 1 ao 15 são 
pontos irradiados, ou seja, os pontos que foram lidos a partir da instalação do equipamento e das visadas. 
Destaca-se que as altitudes dos pontos irradiados são determinadas com base na leitura de ré da estação 
correspondente. 


Convém ressaltar que este método é recomendado para grandes extensões de áreas e quando se observa 
que o terreno possui uma inclinação muito expressiva não sendo viável ter a visualização de todos os 
pontos apenas em uma única instalação do Nível (SILVA; SEGANTINE, 2015). 


No entanto, devido a vários fatores inerentes ao processo de medição, existe uma série de erros que 
podem afetar o resultado do nivelamento. Por esta razão deve-se fechar o caminhamento em um ponto 
com altitude conhecida, podendo ser a mesma ou diferente do ponto de partida. Desta forma, pode-se 
obter um parâmetro de qualidade do nivelamento (GHILANI; WOLF, 2012; KAVANAGH; MASTIN, 2014). 


Para a realização dos procedimentos matemáticos para determinar o erro de fechamento altimétrico, 
deve-se obter a diferença entre a altitude do ponto conhecida e a calculada (Equação 4) ou comparar os 
desníveis determinados com a diferença das altitudes conhecidas com os desníveis obtidos em campo 
(Equação 5) (VEIGA, ZANETTI; FAGGION, 2012). 


ea = o eenisairado Ei Ds or necuã) (8) 


E À Ahcaiculados E e o a) (5) 67, 


Para que o erro de fechamento possa ser distribuído é necessário que ele seja menor que a tolerância do 
projeto. Normalmente os erros do nivelamento são acidentais em função do comprimento da linha de 
visada (SILVA; SEGANTINE, 2015). Por essa razão, a tolerância do erro de fechamento do nivelamento 
geométrico (T) é expressa pela Equação 6 (ABNT, 1994). 
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T = kvd (6) 


Onde: 
k: constante adotada em função da precisão considerada no nivelamento; 


d: comprimento total do trecho nivelado, em quilômetros. 


Segundo Garcia & Piedade (1984), classifica-se em: 


n Alta ordem: o erro médio admitido é de +1,5mm/km percorrido; 

u Primeira ordem: o erro médio admitido é de +2,5mm/km percorrido; 
= Segunda ordem: o erro médio admitido é de 1,0cm/km percorrido; 

u Terceira ordem: o erro médio admitido é de 3,0cm/km percorrido; 

= Quarta ordem: o erro médio admitido é de 10,0cm/km percorrido. 


Nota-se um rigor maior em relação à tolerância ao longo dos anos, especialmente seguindo critérios de 
precisão estabelecidos pela NBR13133 para a execução de atividades desse tipo, especialmente quando se 
trata de obras que envolvem engenharia de estradas, engenharia civil e outras modalidades de engenharia 
que necessitam de resultados muito precisos para a boa execução das mesmas. 


Dentre o que foi exposto, é tido como nivelamento com qualidade aquele que o erro de fechamento (eq) é 
menor que a tolerância (T), caso isso não ocorra deve-se realizar novas observações em campo, 
atendendo-se aos seguintes pontos, que segundo Ghilani e Wolf (2012), são fontes de erro mais comuns no 
nivelamento; 


= Execução sem a devida verificação do nivelamento do nível de bolha do aparelho; 

= Leitura incorretas da mira; 

" Erro de anotação da mira; 

n Erro de verticalidade da mira; 

= Erro da graduação da mira; 

= Erro de horizontalidade do eixo de visada (eixo óptico não está perpendicular com o eixo vertical). 
u Erro devido a influência da curvatura e refração terrestre; 

= Erro na na determinação a altitude dos Pontos de Mudança; 

" Focagem incorreta dos retículos estadimétricos; e 

= Calor ou reverberação e ventos. 


Z 


Por outro lado, caso esteja dentro das tolerâncias pré-estabelecidas, é importante realizar as 
compensações para que os erros se dilua nas observações. A compensação (Ch) é oposta ao erro (ea), ou 
seja, Ch = - ea. Na verdade, Ch é distribuído repartindo o valor a cada ponto, podendo ser empregado 
proporcionalmente ao número de pontos de mudança (Equação 7) ou inversamente proporcional a 
distância do lance (Equação 8) (SILVA; SEGANTINE, 2015). Vala-se ressaltar, que a Equação 8, considera- 
se que quanto maior a distância do lance, maior é o erro. 


C(n=—— (7) 


Cn AC E! (8) 


Onde: 
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n: número de pontos de mudança do caminhamento ou números de desníveis. 
> d: é a comprimento total do trecho nivelado, ou seja, é o somatório das distâncias de cada visada D; 


De modo que, as distâncias de cada visada (Nível à Mira) são determinadas a partir a equação 9, ou seja, é 
a subtração do fio superior e fio inferior, multiplicada pela constante estadimétrica da mira (L), sendo 
normalmente igual a 100. 


Di (nível àmira A) = (FSa — Fla) x L (9) 


Por fim, chega-se a etapa de se realizar a correção das altitudes a partir da compensação. Atém-se que 
nesse método não são definidas coordenadas planimétricas (X, Y) (KAVANAGH; MASTIN, 2014). Logo, a 
partir da Equação (10) determina-se as altitudes corrigidas dos pontos. 


H, =H4 + Ahi;-, duo E (10) 


6.5. PRÁTICA 


Em síntese, os materiais recomendados para esta aula prática serão: 


= 1 Nível Topográfico; 

n 1 Tripé; 

u 2 Mira estadimétrica; 

n 2 Nível de cantoneira; 

u 1 Trena de 30 m; 

= 9 Piquetes; e 

= Prancheta e recursos necessários para anotação de dados. 


6.5.1. PROCEDIMENTOS EM CAMPO 


Recomenda-se adotar grupos de quatro alunos, de modo que dois destes verticalizam as miras 
estadimétricas na Ré e na Vante; outro aluno, portanto, instale, nivele e opere o Nível assim como realiza 
as leituras das miras; e o último auxilia o operador, anotando e verificando as leituras a partir dos cálculos. 
No entanto, vale frisar que essa quantidade é para fins didáticos, pois na prática somente dois 
profissionais são suficientes para esta operação. 


Considerando essa quantidade de alunos por grupo, sugere-se uma poligonal com no mínimo oito pontos 
para um caminhamento fechado (começa e termina no mesmo ponto de referência) ou nove pontos para 
encaminhamento apoiado (começa e termina em pontos de referência diferente), totalizando oito 
estacionadas. Este procedimento é para que os alunos tenham a oportunidade de operar pelo menos duas 
vezes o equipamento, desde que seja feito o devido rodízio de funções entre os operadores. 


O objetivo deste procedimento prático é para que os alunos consigam executar todas as etapas do 
nivelamento geométrico, iniciando com o planejamento dos pontos, partindo com as coletadas de dados, 
realização dos cálculos apresentados no tópico anterior (3.4), bem como analisar o erro e distribuição do 
mesmo. 


Primeiramente, deve-se checar se o ponto de partida e chegada tem altitude conhecida para realizar o 
transporte das altitudes e realizar a materialização dos pontos de interesse com piquetes, chapas ou pinos, 
encravados no solo. Caso julguem necessário, pode-se empregar estacas a fim de auxiliar a localização 
destes pontos. 
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Posteriormente, dois alunos devem colocar as miras sobre os pontos de interesse do primeiro lance. As 
miras deverão estar verticalizadas com o auxílio do nível de cantoneira e serão apoiadas sobre chapas, 
pinos ou piquetes materializados, e no caminhamento, deverá colocar sobre sapatas, evitando apoiar 
diretamente sobre o solo. Vale frisar que neste método, diferentemente dos demais métodos 


planialtimétricos, o Nível não é estacionado nos pontos de interesse, assim, não há necessidade de 
centralizar o Nível em nenhum ponto materializado. 


O Nível deve ser instalado sobre o tripé a meia equidistância das miras, de modo a minimizar o efeito da 
curvatura terrestre e da refração atmosférica, bem como compensar a horizontalidade do eixo de 
colimação da luneta (IBGE, 2017). Para tal fim, recomenda-se empregar uma trena de 30m ou realizar 
uma contagem a passo, já que se aceita uma discrepância entre a distância para a Ré e a distância para a 
Vante de 3m. Além disso, para uma leitura mais precisa recomenda-se que esta distância seja inferior a 
60m (IBGE, 2017). Após a instalação, nivela-o com o auxílio dos calantes e tripé, conforme já visto em sala 
de aula. 


Em seguida, orienta-se um aluno operar o Nível e ler na mira à Ré os fios estadimétricos (FS, FM, FI) e 
outro aluno anotar na planilha (APÊNDICE 1) e verificar a leitura do FM lido está condizente com a 
equação 11, de modo que a diferença entre o FM lido e FM o calculado não deve exceder 2mm. Em seguida 
visar a mira à Vante repetir os procedimentos da ré. Vale ressaltar que as leituras nas miras deverão se 
situar acima de 20-50 cm do solo para evitar turbulências causadas pela reverberação (IBGE, 2017). 


FM+FI 
FMcaiculado = 2 (11) 


Ao realizar o cálculo do desnível do primeiro lance, deve-se seguir para o próximo lance tomando o 
cuidado de alternar as miras, ou seja, a mira que estava posicionada no ponto Ré seja posicionada neste 
lance no ponto Vante. Este processo possibilita minimizar um possível erro de índice da Mira (IBGE, 
2017). 


Este processo deve ser feito até o último lance da poligonal, com alternância das funções dos alunos de 
modo que todos os alunos do grupo passem por todas as etapas (instalar e nivela o Nível, aprumar a mira, 
realizar as leituras e os cálculos) afim de fixar os conhecimentos. 


Após o procedimento de campo e coletar os dados, deve-se finalizar os cálculos dos desníveis (Equação 1), 
do transporte das altitudes (Equação 2), erro de fechamento altimétrico (Equação 4 ou 5), conforme 
descrito no tópico 3.4. 


Caso o erro de fechamento altimétrico seja maior que a tolerância estabelecida (1,0cm/km), o trabalho 
deverá ser refeito, tomando-se os devidos cuidados enfatizados no item anterior (seção 3.4). Caso 
contrário, deverão proceder com os cálculos de compensação do erro altimétrico (Equação 7 e 8) e 
determinação da altitude (Equação 10). 
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APÊNDICE 1 
NIVELAMENTO GEOMÉTRICO 
Saída = ;H= ; Chegada = :H= 


Ponto Leitura Mira Distância Altitude 


Leitura Mira Ré Distância Ré AH | Compens.C 


Visado Vante Vante Corrigida 


Transp. Lam 


RH —— 


E o a 
E A [O O 
Bo To TJmp o DO 


E PR | SE | DS RR EE 
2d= 


Capítulo 7 


Aula prática: Sistema de recalque 
Levantamento da curva de uma 
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7.1 INTRODUÇÃO 


A necessidade de transportar água de um reservatório situado em uma cota inferior para outro situado em 
uma cota superior, impõe a existência de um equipamento eletromecânico que forneça uma certa 
quantidade de energia ao líquido. Esta energia é expressa em termos de coluna de líquido e é denominada 
altura manométrica total (Hman), dada por: 


Hman = Hg + perdas = Hr+Hs (1) 
Onde: 


Hg - é a altura geométrica, diferença de níveis entre os reservatórios. (Hr+Hs) 


Hr - é a altura manométrica de recalque, diferença entre a cota piezométrica na saída da | bomba ea cota 
do eixo. 


Hs - é a altura manométrica de sucção, diferença entre a cota piezométrica na entrada da bomba e a cota 
do eixo. 


A altura manométrica de recalque é sempre positiva e a altura manométrica de sucção pode ser negativa 
ou positiva, dependendo se a linha piezométrica corta ou não a tubulação de sucção. No caso da instalação 
montada a bomba encontra-se não afogada, ou seja, o eixo da bomba está acima do nível do líquido. 


As perdas de carga localizadas, ocasionadas pela mudança de seção do escoamento ou direção do fluido, 
diferenciam-se das perdas distribuídas, uma vez que estas dependem de parâmetros definidos e bem 
conhecidos, enquanto aquelas dependem de fatores de difícil determinação e bastante diversificados. As 
perdas de carga localizadas variam de acordo com as características do escoamento à montante e jusante 
da singularidade e, apesar de denominadas localizadas, na realidade se verificam em um pequeno trecho 
de um lado ou de outro do acessório. Em outras palavras, a queda da linha de energia, devido a presença 
do acessório, não ocorre verticalmente. 


Para tubulações de PVC, de diâmetros menores que 4”, transportando água fria, a perda de carga unitária 
pode ser calculada pela equação de Fair-Whipple-Hsiao: 


1,75 


ne o ioigss SL om mm OU e) 
DEE 


7.2 DESENVOLVIMENTO 
7.2.1 OBJETIVO GERAL 


Levantar experimentalmente a curva de uma bomba centrifuga e comparar com a curva do fabricante. 


7.2.1.1 OBJETIVOS ESPECÍFICOS 
= Operacionalizar um experimento de medição de vazão; 
= Observar o funcionamento de um medidor de vazão; 
= Identificar as partes constituintes de um sistema de turbobombas; 


= Confeccionar a curva da bomba com os dados obtidos; 


7.2.2 METODOLOGIA 


No laboratório de hidráulica do curso de engenharia sanitária da UFPA, está instalado um sistema de 
recalque com duas bombas em paralelo marca KSB, com manovacuômetro na entrada da bomba e 
manômetro (ambos marca Mecaltec) na saída da mesma. O sistema é composto com tubulações de PVC, 
sendo 75 mm na sucção e 50 mm no recalque. A bomba encontrava-se não afogada, e succiona do 
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reservatório 1 de 500 L, para o reservatório 2 de 310 L, sendo que deste a água retorna para o 
reservatório 1 por gravidade, ver Figura 1. 


Figura 1- Instalação de recalque com recirculação 


Reservatório 2 


Reservatório 1 


Na Figura 2 verifica-se detalhes da bomba e dos medidores de pressão. 


Figura 2-Bomba centrífuga e manômetros de medição 


Na Figura 3 indica-se o medidor de vazão (de engrenagens) oval Tech-meter, que foi utilizado na pesquisa, 
e as adaptações feitas para o funcionamento do mesmo na tubulação de recalque do sistema de bombas. 
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Figura 3- Medidor de engrenagens (vazão em L/min) 
RA de iSTi e 


da tá 


Figura 4- Detalhe do medidor na tabulação 


Procedimentos práticos desenvolvidos: 


1. Realizar os levantamentos, a trena, das cotas lineares e conexões para a obtenção dos 
comprimentos reais e equivalentes para determinação da perda de carga; 


2. Elaborar o desenho isométrico da instalação; 
3. Escorvar e ligar a bomba; 
4. Conferir as condições de bombeamento; 


5. Após estabilização da vazão, para dez diferentes manobras de registro na caixa dágua de 310L, 
anotar os dados de vazão para 10 medições, no medidor de engrenagens (Flowpet - OVAL 
TechMeter); 


6. Para cada medida de vazão anotar as pressões correspondentes na saída e entrada da bomba; 
7. Construir o gráfico do ensaio (curva da bomba); 
8. Comparar com dados do fabricante. 


Pela tabela abaixo quantificou-se os comprimentos equivalentes das conexões e determinou-se as 
perdas de carga. 
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Tabelal- Comprimentos equivalentes 


Perdas de carga localizadas - Sua equivalência em metros de tubulação de PVC rigido 


Joelho Joelho Curva Curva Tê90* Té90* TE90* Entrada Entrada Saídade Válvula de Válvula de Válvula de 


se as 9” 45º Passagem Saída Saída Normal de Borda Canalização PéeCrivo Retenção Retenção 
Direita delado Bilateral Tipo 


Co voto re 


Fonte: Tigre manual técnico 


7.2.3 RESULTADOS DE UM EXPERIMENTO: 


Os resultados das medições efetuadas e dimensionamentos realizados são apresentados nas tabelas 
abaixo. 
Tabela 2- resultados dos 10 experimentos. 
Vazão Vazão Vacuômetro Manômetro Hman 
(L/min) (m?/s) (mca) (mca) (mca) 


Ensaio 


Para tubulações de PVC, de diâmetros menores que 4”, transportando água fria, a perda de carga unitária 
pode ser calculada pela equação de Fair-Whipple-Hsiao: 


J = 0,00086 ———= com J(m/m), O(m [e D(m) 
p4:15 
2) B]=[2).B) 
2=PRR+A RR 


Deste modo calculou-se as perdas de carga para cada vazão verificada, na sucção e no recalque do sistema. 
Pelo levantamento in loco foi determinada comprimento real da sucção de 0,85 m. O comprimento 
equivalente foi referente a uma curva 90º (1,4 m) e uma válvula de pé com crivo (25 m). (Tabela 3) 
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Tabela 3- Perda de carga na sucção D=75mm 
Ds (m) Q (m3/s) | J Lreal Leg Hfs(m) 


No recalque mediu-se 11, 31 m de comprimento linear e 39,1 m de comprimento equivalente, referentes a 
2 válvulas de gaveta, 1 Tê saída lateral, 1 tê passagem direta, 7 joelhos 90º, 02 joelhos 45º. 


Tabela 4-Perda de carga no recalque D=50 mm 


Dr (m) Q (m3/s) Lreal Leq Hfs(m) 


Na tabela 5 apresenta-se a perda total do sistema. 


Tabela 5-Perda total no sistema. 
Perda 


Ensaio Vazão (L/min) Vazão(m?/s) Hman(m) Hfr(m) Hfs(m) total (m) 


1] Sa | omoMos | das | Om | cora | 035 | 
2 | sa | omoms | no | 02 | aos | amu | 
[= [55 [omnes | us | ozma | aos | ozom | 
[= | 6 | omom | 124 | 025 | ao | az | 
[5 [4 [omorimasa | sa2 | oz | 005 | oz | 
6 | ser [omni | 1 | om | oo | axis | 
[e | ssa [om | 156 | 01 | ooo | orcs | 
[5 [227 [omosamaa | ama | 0105 | ooma | anã | 


Nas Figuras 5 e 6 são representados, respectivamente, a vista frontal do sistema proposto e a curva 
característica da bomba levantada experimentalmente. 
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Figura 5- Detalhes do sistema 


VISTA FRUNTAL 


Figura 6- Curva característica experimental da bomba. 
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Comparando-se as Figuras 6 e 7, percebe-se que o rotor que mais se aproxima da curva da bomba 
levantada experimentalmente foi o de 195 mm, para o modelo 25-200 com 1750 rpm. 
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Figura 7- curva do fabricante. 
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Capítulo 8 


Aula prática: Carga fixa-poluição 


Neyson Martins Mendonça 
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8.1. INTRODUÇÃO: QUALIDADE DA ÁGUA E OS TIPOS DE FONTES 


A água é usada para diversos fins, como consumo humano, lazer, irrigação, entre outros e nos dias atuais é 
considerada como uma “commodity” preciosa e cuja escassez tem sido tema de preocupação em diversos 
países em razão da disponibilidade dessa em termos de quantidade e de qualidade, frente as alterações 
climáticas e situações de desastres e emergências que possam acontecer. 


A ANA menciona que no Brasil o conhecimento sobre a qualidade das águas está entre um dos mais 
relevantes e emblemáticos desafios a serem trabalhados com a devida atenção na área ambiental desse 
Pais, em decorrência de se ter informações esparsas ou inexistentes quanto ao monitoramento da 
qualidade dos corpos hídricos brasileiros, capazes de subsidiar um sistema de gestão hídrica fundamento 
numa economia ambientalmente sustentável e socialmente justa. 


Um dos principais objetivos para implantação desse sistema de gestão hídrica permeia pela identificação 
das fontes pontuais que são locais fixos de descarga dos efluentes e das fontes não pontuais, também 
conhecidas como difusas que podem ocorrer devido ao escoamento superficial ou pela transferência da 
atmosfera para o sistema aquático. 


Seja para fontes pontuais ou não pontuais, o que se deseja caracterizar é a variabilidade espaço-temporal 
da carga liberada para o corpo d'água receptor, a qual é normalmente expressa em unidades de massa por 
unidade de tempo para um dado conjunto de poluentes químicos, físicos e biológicos. 


Assim no presente capítulo será tratado o procedimento técnico prático para se obter a carga de poluição 
oriunda de fonte pontual de descarga de efluente, tendo por base as notas de aulas da Disciplina TE06144 
- Qualidade da Agua do Curso de Engenharia Sanitária e Ambiental da UFPA. 


8.2. OBJETIVOS INSTRUCIONAIS 


Ao final desta aula prática de caracterização de fonte fixa de poluição para descarga dos efluentes o 
discente deverá ser capaz de: 


1). Identificar os locais para medição da descarga dos efluentes (bruto ou tratado); 
2). Identificar os principais dispostos e/ou equipamentos para caracterização de fonte fixa de poluição; 


3). Obtenção do volume representativo de amostra composta para caracterização de fonte fixa de 
poluição; 


4). Organizar e tratar os dados mensurados em campo para se obter a descarga de poluição para fonte 
pontual fixa. 


8.3. MÉTODO DE MEDIÇÃO DA CARGA FIXA DE POLUIÇÃO 


O método de medição da carga fixa de poluição segue as diretrizes da NBR 13402 da ABNT (1995) que 
versa sobre caracterização de cargas poluidoras em efluentes líquidos industriais e domésticos, e que para 
aplicação desse necessita que se observe sobre os seguintes aspectos principais, a saber: 


-Demarcação de ponto de coleta: O ponto de coleta deve estar situado necessariamente no efluente da 
fonte geradora da carga poluidora; 


-Determinação de vazão: deve ser efetuada utilizando-se dispositivos adequados, como: calha Parshall, 
calha palmer bowlus, e medidor de vazão eletromagnético. 
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Tabela 1. Principais equipamentos de medição de vazão de acordo com o tipo de escoamento em conduto 
forçado ou livre. 


Equipamento de medição de vazão 


Desenho técnico 


Medidor de vazão eletromagnético (MQE): esse equipamento aplica o princípio de 
medição de vazão com base na lei de Faraday que indica a tensão do sinal (E) é 
dependente da velocidade média do líquido (v) da intensidade do campo magnético 
(B) e do comprimento do condutor (D) associada a equação da continuidade Q = A. v, 
para mensurar a vazão em m3/h tubulações; 


Calha Palmer Bowlus: esse dispositivo de medição de vazão é utilizado para tubos 
parcialmente cheios, o qual pode ser acoplado diretamente na tubulação de líquidos 
fluindo por gravidade para em conjunto com medidor de nível ultrassônico registrar a 
distância entre o transdutor e a superfície monitorada (lâmina de água), cujo valor é 
posteriormente inserido em formulação para obtenção da vazão (m3/h); 


Calha Parshall: esse dispositivo consiste numa seção convergente, numa seção 
estrangulada “garganta” e uma seção divergente, dispostas em planta. O escoamento 
do fluido submetido a essa unidade é tranquilizado em sua seção convergente 
evitando os efeitos da velocidade, a partir de então com auxílio de medidor de nível 
ultrassônico se registra a distância entre o transdutor e a superfície monitorada 
(lâmina de água), cujo valor é posteriormente inserido em formulação para obtenção 
da vazão (m3/h). 


-Definição das variáveis de poluição de interesse a serem determinadas: considerando a fonte pontual 
quanto a origem doméstica ou industrial e a legislação ambiental brasileira em termos do conjunto de 
poluentes químicos, físicos e biológicos para se caracterizar a carga poluidora gerada, recomenda-se a 
princípio realizar a análise de todos os contaminantes prováveis de serem encontrados, mediante consulta 
a Tabela 8 do Guia nacional de coleta e preservação de amostras CETESB (2011) e/ou a Tabela 1 da NBR 
13402 da ABNT (1995). A 


Tabela 2 a seguir ilustra exemplo de agrupamento de variáveis físico-químicas e os exames 
microbiológicos para monitoramento de lixiviado em aterro sanitário 


Tabela 2. Grupo de variáveis físico-químicas e os exames microbiológicos para monitoramento de lixiviado 
em aterro sanitário. 


Grupo de variáveis físico-químicas e os exames microbiológicos Legislação e/ou Referência técnica 


Resolução CONAMA Nº 430/2011 
(Art 16. Condições de Lançamento) 


Grupo 0: pH, Temperatura, SD (sólidos sedimentáveis), 
DBOs,oc, OG (minerais) e OG (vegetais) 


Grupo 1: Arsênio (As), Bário (Ba), Boro (B), Cádmio (Cd), 
Chumbo (Pb) 

Cobre (Cu) dissolvido, Estanho total (Sn), Fe (Fe) dissolvido, Mn 
(Mn) dissolvido, Mercúrio (Hg), Prata (Ag), Selênio (Se), Níquel 
(Ni), Nitrogênio amoniacal total, Prata total (Ag), Selênio total 
(Se), Sulfeto (S) e Zinco (Zn) 


Resolução CONAMA Nº 430/2011 
(Tabela I- Condições de Lançamento -Parâmetro 
inorgânico) 


Grupo 2: Benzeno, Clorofórmio, Dicloroeteno, Estireno, Resolução CONAMA Nº 430/2011 


Etilbenzeno, Fenóis totais, Tetracloreto de carbono, 
Tricloroeteno, Tolueno e Xileno 


(Tabela I- Condições de Lançamento -Parâmetro 
orgânico) 


Grupo 3: ST (Fixo e volátil), SST (Fixo e volátil), Condutividade, 
AT, Cor aparente, Turbidez, Dureza Total (Ca2+; Mg2+), NTK, P, 
Sulfato, Cloreto (Cl-), DQOs, DQOs Cloro* e THM* 
Nota:*Quando se faz cloração no efluente tratado 


Usuais para controle da ETE e tem valores de 
referência estabelecidos para níveis de restrições 
dos padrões de lançamento (DAE artigo edição 215 
Nº 1764) 


Grupo 4: Coliformes totais e E.coli; 
Nota: Em relação ao Grupo 4, critérios relativos a reuso 
também irão ser observados. 


Resolução CONAMA Nº 357/2005 
Resolução CONAMA Nº 274/2001 


- Finalizado a definição das variáveis de poluição de interesse a determinação da carga fixa de poluição, 
deve ser realiza para cada variável, utilizando-se o equacionamento a seguir: 


Per=C.Q.f 


(01) 


"A 
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Em que: 

Pcr: carga poluidora fixa, em kg/dia; 
€: concentração, em mg/L; 

Q = vazão, emm3/s 


f = fator de conversão de unidade (86,4) 


De acordo com a NBR 13402 da ABNT (1995) quando se trabalhar com um número (n) reduzido de 
amostra (mínimo de seis para cada ponto), a determinação da carga poluidora Pcr deve ser calculada 
considerando as seguintes equações: 


a (Q + —) (c + E) (02) 

im (q E —e) : [ E o) (03) 

use [O 
Em que: 


Pcr: carga poluidora fixa; 

Q: média aritmética de vazão; 

So: desvio padrão para vazão; 

C: média aritmética de concentração; 

Sc: desvio padrão para concentração; 

t: coeficiente de Student (ver Tabela 2 da NBR 13402 da ABNT (1995); 


n: número de amostra. 


- E para a obtenção do volume representativo de amostra para caracterização de fonte fixa de poluição, se 
utiliza o equacionamento previsto pela CETESB (2011), conforme se pode observar a seguir: 


Vais ( = ) (05) 


Em que: 
Q;: vazão do horário (L/s); 
Q: vazão média (L/s); 
n=numero de determinações; 
Var=volume de alíquota de amostra correspondente ao horário (mL); 


Vac=volume de amostra composta (mL). 


Práticas em Saneamento, Clima e Meio Ambiente 


8.4. ANTES DA MEDIÇÃO DA CARGA FIXA DE POLUIÇÃO 


Deve-se separar os seguintes materiais antes de se deslocar para o local de medição: 


21) Protetor solar; 

22) Tenda ou guarda sol e cadeiras; 

23) Planilha de campo (Apêndice 1), prancheta, caneta e lapiseira; 

24) Identificação no local de medição do dispositivo de vazão que será utilizado para media da 


descarga: calha Parshall, calha palmer bowlus, e medidor de vazão eletromagnético; 


25) 50 Frascos PP (plástico de polipropileno autoclávavel) com graduação com tampa rosqueável 
antivazamento para armazenamento de 250mL devidamente rotulados e identificados; 


26) Caixa isotérmica com dimensional de 39,5 cm x 27,5 cm x 34 cm e gelo; 
27) Termômetro Infravermelho Digital de LCD Portátil para faixa de medição de -2ºC a 40ºC; 


28) Par luva de nitrilo preta sem pó descartável (Tamanho: P,M ou G); 


Não esqueça de estar vestido adequadamente de calça, botas, camisa confortável, óculos de proteção solar 
e chapéu/boné. Também fiquem atentos para fazer registros fotográficos em todas as etapas. 


8.5. ATIVIDADES EM CAMPO 


Nesta atividade de campo será mensurada a descarga do efluente bruto e o volume representativo de 
amostra composta para caracterização de fonte fixa de poluição, para isso inicialmente serão formados 04 
grupos contendo cinco 5 alunos para numa dada estação de tratamento de esgoto (ETE), identificarem o 
local de chegada da descarga do efluente e o equipamento de medição de vazão utilizado nessa unidade, 
mediante a elaboração do fluxograma das etapas que ocorrem nessa unidade tanto para a fase líquida 
como sólida. 


A seguir, os grupos de alunos serão divididos por turnos de medição, obedecendo a seguinte ordem: Grupo 
1-7:00 h até 9:30 h; Grupo 2-10:00 h até 12:30 h; Grupo 3-13:00 h até 15:30 h e Grupo 4-16:00 h até 18:30 
h. 


Cada grupo de medição, deverá em seu turno, realizar seis (06) leituras de medição de vazão e a obtenção 
mediante coleta simples de seis (06) amostras do efluente bruto com auxílio de coletor manual, para 
transferência dessas para os frascos PP de 250 mL, acondicionados com gelo na caixa isotérmica, para 
posterior envio das amostras para laboratório de análises e determinação da demanda química de 
oxigênio (DQO). 


Em cada grupo de medição, um (01) aluno deverá ser responsável pelo registro de fotográfico e 
preenchimento da ficha de campo, dois (02) alunos deveram ficar dedicados a medição de vazão e os dois 
(02) alunos restantes ficam envolvidos com a coleta de amostras do efluente bruto para obtenção da DQO. 


8.6. TRATAMENTO DOS DADOS OBTIDOS EM CAMPO 


Todos os dados obtidos nas atividades em campo devem ser organizados em planilha eletrônica Excel ou 
equivalente e trabalhados para a determinação da carga fixa de poluição e do volume representativo de 
amostra composta para caracterização de fonte fixa de poluição, mediante o uso do equacionamento 
descrito no item 8.3 do capitulo em questão. 


E utilizando o Microsoft Excel ou Microsoft Minitab o tratamento de dados deverá conter os seguintes 
itens principais: 


a) Estatística descritiva para vazão e DQO: número de determinações (n), máximo (máx), médio (x), 
mínimo (mín), desvio padrão (DP) e coeficiente de variação (x/DP); 


b) Construir gráficos do comportamento temporal da vazão e da DQO; 


c) Determinação da carga poluidora fixa (média e mediana), máxima e mínima; 
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d) Obtenção do volume representativo de amostra composta, considerado as seguintes condições: Grupo 
1-1.500 mL; Grupo 2-2.730 mL, Grupo 3-3.850 mL.e Grupo 4-835 mL. 


e) Construir gráfico do comportamento temporal do volume de alíquota de amostra correspondente ao 
horário; 
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APÊNDICE 


UNIVERSIDADE FEDERAL DO PARÁ 
INSTITUTO DE TECNOLOGIA- FACULDADE DE ENGENHARIA SANITÁRIA E AMBIENTAL 
LABORATÓRIO DE INSTALAÇÕES PILOTO DE TRATABILIDADE DE ÁGUA E DE LODO 
DISCIPLINA TE06144 - QUALIDADE DA ÁGUA 
FORMULÁRIO DE RELAÇÃO DE EQUIPAMENTO E MATERIAIS PARA TRABALHO DE CAMPO 


Instituição onde será o trabalho de campo: 


LOCALIZAÇÃO GEOGRÁFICA: Instituição responsável pelo material e equipamento: 


LATITUTE: á o há LONGITUTE: *| Contato: Celular: 


E-mail: 


IDENTIFICAÇÃO DE MATERIAIS E EQUIPAMENTOS DE TRABALHO DE CAMPO 


Descrição de materiais e/ou equipamentos Descrição de materiais e/ou equipamentos 


Rua Augusto Corrêa 01 Cidade Universitária CEP 66075-110 Guamá - Caixa postal 479 PABX +55 91 3201-7000 — Belém- Pará Brasil - Universidade Federal do Pará- 


UNIVERSIDADE FEDERAL DO PARÁ 
INSTITUTO DE TECNOLOGIA- FACULDADE DE ENGENHARIA SANITÁRIA E AMBIENTAL 
LABORATÓRIO DE INSTALAÇÕES PILOTO DE TRATABILIDADE DE ÁGUA E DE LODO ITEC;: 
DISCIPLINA TE06144 - QUALIDADE DA ÁGUA 
FORMULÁRIO DE AMOSTRAGEM DE CAMPO DE COLETA DE AMOSTRA DE ÁGUA E DE LODO 


TIPO DE COLETA: 


LOCALIZAÇÃO GEOGRÁFICA: INSTITUIÇÃO RESPONSÁVEL PELA COLETA: 


LATITUTE: hz LONGITUTE: ASSINATURA: 


CONDIÇÕES DO TEMPO: 


IDENTIFICAÇÃO DA AMOSTRA E DESCRIÇÃO DO LOCAL DE AMOSTRAGEM 


VOLUME DE TIPO DE 


ID DA AMOSTRA DESCIRÇÃO DO LOCAL DE AMOSTRAGEM TEMPO 
AMOSTRA (mL) | AMOSTRA 


TIPO DE AMOSTRA: AC = água de consumo, AT = água tratada, AB-=água bruta; EB= esgoto bruto; ET=efluente tratado; LAN=lodo anaeróbio; 
LAE=lodo aeróbio; El=efluente industrial 


OBS: 


Rua Augusto Corrêa 01 Cidade Universitária CEP 66075-110 Guamá - Caixa postal 479 PABX +55 91 3201-7000 — Belém- Pará Brasil - Universidade Federal do Pará- 
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UNIVERSIDADE FEDERAL DO PARÁ 
INSTITUTO DE TECNOLOGIA- FACULDADE DE ENGENHARIA SANITÁRIA E AMBIENTAL 
LABORATÓRIO DE INSTALAÇÕES PILOTO DE TRATABILIDADE DE ÁGUA E DE LODO ITE C Diria 
DISCIPLINA TE06144 - QUALIDADE DA ÁGUA 
EXEMPLO DE FLUXOGRAMA DAS ETAPAS DE TRATAMENTO DA ETE-FASE LÍQUIDA E SÓLIDA - PARA MAPEAMENTO DE CARGA 
- DESINFECÇÃO - “TRATAMENTO SECUNDÁRIO- 
METANO 


(Queima + secagem), — —p— 


-Disposição Final. 


-Acondicionamento- 
em container de PEAD Sólidos [Transporte ) 
com cal (Ca(OH)2) DR Ee ga —+> 
tempo recomendado 
o ATERRO SANITÁRIO 
4 


di 


z 


Oo 


Esgoto Desarenado 


Prniciar = 62.602 hab 
Prna =73.927 hab 


EUA do 


DESINFECTADO 


FINAL Quicia = 115 Us 


Qua =131 Lis 


EN 


Extensão da rede = 54.374m 


-Acondicionamento- 
NÃo em container de PEAD 
> comcal(Ca(0H)2) — — D [Transporte 
excesso E tempo recomendado “+ UMa + ] 
- TRATAMENTO DO LODO - I 
(Condicionamento + Desaguamento + Higienização) Tratamento- -Disposição Final- ] 
] Areia deverá sofrer | 
I 
I 
| 


Uso para produção de 
1 4 
osicão Find — — — > "ente — — —p agregado no setor da 
-Disposição Final. peneirada construção civil 


-Projeto de paisagismo de área verde do município de Belém (PA) 
> | -Jardinagem de empreendimentos habitacionais de interesse I 
lodo seco y social 
Á Condicionador de solo 
- agricultura familiar na área insular, 
—> plantação de grama 
- recuperação de área degradada 
- revegetação de ADRS, etc. 


Rua Augusto Corrêa 01 Cidade Universitária CEP 66075-110 Guamá - Caixa postal 479 PABX +55 91 3201-7000 — Belém- Pará Brasil - Universidade Federal do Pará- 
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INSTITUTO DE TECNOLOGIA 
FACULDADE DE ENGENHARIA SANITÁRIA E AMBIENTAL ITE C: UEra” 
LABORATÓRIO DE INSTALAÇÕES PILOTO DE TRATABILIDADE DE ÁGUAS E LODO 


DISCIPLINA TE06144 - QUALIDADE DA ÁGUA 
Unidade: DATA 
Local de medição: 
Método de medição da vazão: 
Tipo de amostragem: composta 


Assinatura do responsável pela medição e determinação: 


UNIVERSIDADE FEDERAL DO PARÍ RS 


DQO, | DQO; DQoO 
rt oo | ct [a | qo | o Lo tr [ | tur 


Preencher os campos em cinza para se obter as informações desejadas 


Rua Augusto Corrda 01 Cidade Universitária CEP 88075-110 Guamá - Calxa poctal 478 PABX +65 81 3201-7000 — Belém- Pará Bracll - Universidade Federal do Pará 
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UNIVERSIDADE FEDERAL DO PARÁ 
INSTITUTO DE TECNOLOGIA AS 
FACULDADE DE ENGENHARIA SANITÁRIA E AMBIENTAL ITEC Fo 
LABORATÓRIO DE INSTALAÇÕES PILOTO DE TRATABILIDADE DE ÁGUA E DE LODO 
DISCIPLINA TE06144 - QUALIDADE DA ÁGUA 


IDENTIFICAÇÃO E REGISTRO DO LOCAL DE AMOSTRAGEM PARA AULA PRÁTICA: CARGA FIXA-POLUIÇÃO 


REGISTRO DO LOCAL DE AMOSTRAGEM PARA LEVANTAMENTO DA CARGA FIXA DE POLUIÇÃO 


E SE 
Coordenada georeferenciada 


Descrição do local de mediação e amostragem: 


Rua Augusto Corrêa 01 Cidade Universitária CEP 66075-110 Guamá - Caixa postal 478 PABX +55 91 3201-7000 — 
Beléêm- Pará Brasil - Universidade Federal do Pará- 2008 


Capítulo 9 


Aula prática: Ensaio de adensamento 
por gravidade em batelada de lodo de 
eta 


Luiza Girard Teixeira 
Marcia Uchôa Matos 
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9.1. OBJETIVOS INSTRUCIONAIS 


Ao final desta aula prática de ensaio de adensamento por gravidade em batelada de lodo de Estação de 
Tratamento de Agua - ETA você deverá ser capaz de: 


1) Compreender o conceito de adensamento por gravidade de lodo de ETA e as variáveis envolvidas no 
processo; 


2) Realizar o ensaio de adensamento por gravidade em batelada de lodo de ETA; 


3) Avaliar os resultados obtidos. 


9.2. ADENSAMENTO POR GRAVIDADE DE LODO DE ETA 


Segundo a NBR 10004 (ABNT, 1987) os lodos gerados nas ETAs são classificados como resíduos sólidos e, 
portanto, devem ser disposto sem ocasionar danos ao meio ambiente. 


As características dos lodos de ETAs podem variar com o processo de tratamento e depende de fatores 
como: qualidade da água bruta, tecnologia de tratamento, características da coagulação, floculação e 
filtração, uso de oxidante, método de limpeza de decantadores e filtros, entre outros (DI BERNARDO et al., 
2011). 


O tratamento desse lodo visa justamente obter um estado sólido ou semi-sólido, removendo a água para 
concentrar os sólidos, diminuindo seu volume e assim gerar condições adequadas para a sua disposição 
final (RICHTER, 2001). A primeira etapa para o tratamento deste resíduo é denominada de adensamento, 
que consiste na concentração de sólidos presentes no lodo, ou seja, visa remover o máximo de água 
possível do lodo, preparando-o para a etapa posterior de desidratação. 


O adensamento de lodos de ETAs pode ser executado por gravidade, por flotação com ar dissolvido ou por 
meio de adensadores mecânicos do tipo centrífugas ou esteira. Cada um apresenta vantagens e 
desvantagens e devem ser levados em consideração seus aspectos econômicos e técnicos para a 
determinação de parâmetros de projeto de cada processo unitário. Assim, a concepção e o 
dimensionamento corretos do sistema de adensamento são de suma importância no sucesso operacional 
do sistema de tratamento do lodo. 


O adensamento por gravidade é o mais utilizado devido a simplicidade de operação e ao baixo consumo de 
energia. Além disso, pode servir como uma unidade de equalização. 


Os adensadores podem ser operados de forma contínua ou por batelada. São normalmente unidades 
circulares com um mecanismo de raspagem de lodo no fundo. Em adensadores contínuos, o lodo entra 
pela parte central e se distribui radialmente. O sobrenadante é retirado pela parte superior através de 
vertedores. O fundo deve possuir uma inclinação de modo a facilitar a coleta de lodo adensado. Os 
adensadores por batelada são frequentemente equipados com tubulações de descarga e normalmente dois 
tanques operam em paralelo. O lodo é colocado em uma unidade até atingir o nível de operação. Após um 
determinado período, o sobrenadante é retirado, e o lodo é bombeado para posterior tratamento e 
disposição final. O sobrenadante do processo de adensamento pode ser retornado ao início da ETA. 


Os lodos gerados nas ETA's são, em sua maioria, de difícil desidratação, tornando-se necessário o seu 
condicionamento para a melhoria de suas características de adensamento e posterior desidratação. O 
condicionamento químico é o mais utilizado, tanto para o processo de adensamento quanto para o de 
desidratação. O processo pode ocorrer pela adição de cloreto férrico, cal ou polímeros. O tipo e a dosagem 
de condicionante variam com a qualidade da água bruta, com o tipo de tratamento e com o processo de 
desidratação a ser utilizado (AWWA, 1996). 


Os polímeros são os produtos mais utilizados no condicionamento de lodos provenientes de ETA's. A 
popularidade dos polímeros se deve à sua facilidade de manuseio, necessidade de pequeno espaço para 
armazenamento e sua eficiência. Podem ser catiônicos, não iônicos ou aniônicos, de acordo com a carga de 
seus grupos funcionais, e de baixo, médio ou alto peso molecular. Os mais utilizados para o 
condicionamento são os de alto peso molecular, e atuam pelo processo de formação de pontes 
interparticulares (AWWA,1989). YUZHU (1996) afirma que a carga tem a mesma importância do peso 
molecular, e que a seleção do polímero depende de testes com o lodo específico. Na Tabela 1 podem-se 
observar valores típicos de dosagem de polímeros. 
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Tabela 1 - Dosagens típicas de polímeros usadas no condicionamento 
Tipo de lodo Dosagem (g/kg ST) Referência 


Lodos de Ale Fe 1-10 Montgomery (1985) 
Lodos de Al 1-10 Yuzhu (1996), Matos et al. (2013) 
Lodos de Fe 2-8 Teixeira (1999), AWWA (1996) 


Para dimensionamento dos adensadores por gravidade necessita-se definir a Velocidade de Sedimentação 
e Zona (VSZ), ou seja, a velocidade da interface na primeira fase da sedimentação e das condições do 
adensamento até o teor de sólidos totais desejados (DIBERNARDO et al. 2002). As relações mais usadas 
para se prever o comportamento da VSZ em relação à concentração de sólidos são: 


V=Ae* Eq.(1) 
W= Axo Eq(2) 


onde: A e B= constantes; 

x = Concentração de sólidos (M.L); 

V = velocidade de sedimentação em zona (L.T-"); 
e = base neperiana. 


A equação (1) foi citada por GREGORY (1979) apud AWWA (1990) e a equação (2) por MONTGOMERY 
(1985). De acordo com a AWWA (1990), a equação (1) é a mais usada e a que tem apresentado melhores 
resultados. 


A obtenção dos coeficientes A e B é feita através de ensaios de sedimentação por batelada, objeto desta 
aula prática. A coluna de sedimentação deve ter pelo menos im de altura e 10 cm de diâmetro, de modo a 
reduzir a discrepância entre os valores encontrados em laboratório e os em escala real. 


A coluna após preenchida com lodo é deixada em repouso, e são registrados os valores da altura da 
interface ao longo do tempo. O procedimento é repetido para os teores de sólidos. A velocidade de 
sedimentação deve ser calculada como a tangente do trecho inicial retilíneo, que terá valor negativo, já que 
o movimento da interface é descendente. 


De posse desses valores obtêm-se um gráfico da velocidade de sedimentação em função da concentração 
de sólidos. 


A equação da velocidade de sedimentação é calculada através da linearização da equação (1) resultando 
na equação (3): 


InV=InA-Bx Eq.(3) 


Os coeficientes A e B são determinados através da curva de regressão linear dos pontos obtidos nos 
ensaios de sedimentação em batelada. 


9.3. ENSAIO DE ADENSAMENTO POR GRAVIDADE EM BATELADA DE LODO DE ETA 


Nossa aula prática acontecerá no laboratório Multiusuário de Tratabilidade de Águas -LAMAG e no 
Laboratório de Instalações Piloto- LAPIL, coordenado pelo Grupo de Estudos em Gerenciamento de Águas 
e Reúso de Efluentes - GESA, localizados no Laboratório de Engenharia Sanitária e Ambiental - LAESA da 
Universidade Federal do Pará - UFPA. 


9.3.1. CARACTERIZAÇÃO DO LODO DE ETA COLETADO 


Uma bombona de lodo de ETA com 5 Litros estará previamente à disposição para a realização da 
caracterização. Para tal, deverá ser realizado o seguinte procedimento: 


Práticas em Saneamento, Clima e Meio Ambiente 


n Homogeneizar o conteúdo das bombonas e retirar 3 amostras de lodo; 


= Determinar o teor de sólidos totais de cada amostra de acordo com o método descrito em APHA, 
AWWA, WEF (2005). 


a Calcular a média aritmética dos três valores obtidos que será designado TST (teor de sólidos 
totais do lodo) 


Nesta etapa você vai precisar: 
3 cápsulas de porcelana de 100 mL, 1 dessecador completo, 1 estufa, 1 proveta de vidro de 100 mL, 1 
pisceta com água destilada e 1 balança de precisão com 4 casas decimais. 


9.3.2. CONDICIONAMENTO DO LODO 


A seleção do polímero será feita seguindo sugestões do fabricante (de acordo com a disponibilidade do 
produto) para uso em lodo de ETA. Assim sendo, será testado um polímero em duas dosagens distintas 
pré-determiandas: 1 e 4 g/kg, que são, segundo a literatura, as dosagens comumente utilizadas no pré- 
condicionamento de lodos de ETAs (TEIXEIRA et al, 1997). O preparo da solução de polímero estará 
descrito no item 1.3.3. 


9.3.3. ENSAIO DE ADENSAMENTO EM BATELADA 


Os ensaios serão conduzidos em coluna de adensamento com 1,00 m de altura e 20 cm de lado (Figura 1). 
A coluna é graduada com a utilização de uma fita métrica fixada à superfície externa da mesma. O volume 
útil da coluna é de 20 litros, fazendo com que a altura do nível d'água seja de 51 cm. Para este trabalho, 
adotou-se um valor médio de G.T (Gradiente de velocidade e tempo de mistura) de 3000, baseado em 
Teixeira (1999), o que resulta na rotação de 115 rpm, de modo a se evitar a quebra dos flocos previamente 
formados durante o pré-condicionamento. Na Figura 2 são mostradas as dimensões da palheta do agitador 
utilizado no ensaio. Serão realizados 3 ensaios de sedimentação de acordo com o estabelecido na Tabela 2. 


Figura 1 - Coluna para ensaio de sedimentação. 
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Figura 2 - Dimensões das palhetas do agitador utilizados durante os ensaios de adensamento 


m=0,/5cm 
rr=3,/5 cm 


Tabela 2 - Resumo dos ensaios de sedimentação em coluna. 


Teor de sólidos na coluna | Dosagem de polímero (g/kg) 
Sem polímero (0,3%) - - 
Com polímero (0,3%) 1 4 


Segundo a metodologia proposta por Teixeira (1999) o ensaio de adensamento em coluna ocorre em 3 
passos, sendo a primeira o preparo da solução de polímero, seguida do preparo do lodo na coluna e por 
ultimo o ensaio de adensamento propriamente dito, conforme apresentado a seguir. 


n Preparo da solução de polímero a 0,5%. 
- Em um béquer com capacidade de 1000 ml adicionar 500 ml de água, medida em balão volumétrico. 


- Com o auxílio de uma seringa, tomar o volume equivalente a 2,5 g de polímero (estado líquido), calculado 
através da equação 1. No caso do polímero em estado sólido pesar diretamente o conteúdo em balança de 
precisão. 


m 
RA Eq. (1) 
pp 

onde: V = volume de polímero (L?) 


m =massa de polímero (M) 


p = massa específica do polímero (M/ Lº) 


- Misturar o conteúdo do béquer com um agitador manual e inserir o volume de polímero previamente 
calculado próximo à região central, de tal forma a não permitir que o polímero fique aderido às paredes do 
béquer ou as hastes do agitador. 


- Homogeneizar durante 60 segundos. 
- Após 15 minutos agitar novamente por mais 60 segundos. 


- Reservar em um balão de 1000 mL. 
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Nesta etapa você vai precisar: 
1 béguer de vidro de 100 mL, 1 balão volumétrico de vidro de 500 mL, 1 balão volumétrico de vidro de 
1000 mL, 1 bastão de vidro, 1 pisceta com água destilada, seringas plásticas de 10 mL, agitador manual 
(tipo mix de cozinha), 1 cronômetro e 1 balança de precisão com 4 casas decimais. 
= Preparo do lodo na coluna 


- Com os valores já determinados de teor de sólidos do lodo coletado, calcular o volume necessário para a 
preparação do lodo, de acordo com o teor de sólidos desejado na coluna, através da equação: 


CMAS Cave Faso) 


onde: Cb = teor de sólidos na bombona (%) 
Vb= volume de lodo a ser disposto na coluna (It) 
Cc= teor de sólidos desejado na coluna (%) 


Vc= volume útil da coluna (It) 


- Despejar o volume de lodo calculado na coluna e preencher o restante do volume com água; 
- Fixar o agitador na parte superior da coluna e homogeneizar o lodo; 


- Retirar uma amostra do lodo da coluna para determinar o teor de sólidos totais de cada amostra de 
acordo com o método descrito em APHA, AWWA, WEF (2005). 


Nesta etapa você vai precisar: 
1 coluna de adensamento (já descrita neste item), 1 cápsula de porcelana de 100 mL, 1 dessecador 
completo, 1 estufa, 1 balde de plástico graduado com capacidade pra 5 L, 1 proveta de vidro de 100 mL, 1 
pisceta com água destilada e 1 balança de precisão com 4 casas decimais. 


n Ensaio de adensamento 
- Ajustar a rotação do agitador (115 rpm), aferindo com um rotâmetro. 


- Inserir, com o auxílio de uma seringa, o volume de polímero calculado de acordo com a dosagem 
desejada. A massa de polímero (mp) deverá ser calculada multiplicando-se a massa de sólidos (kg) 
presente na coluna pela dosagem de polímero (g/kg). 


- Imediatamente após a adição do polímero, deixar em mistura rápida por um tempo de um minuto. 


-Parar a agitação, ligar o cronômetro e anotar a posição da interface na coluna a medida que esta desce, em 
intervalos de tempos regulares de 30 segundos, até que tenham decorrido dez minutos. 


- Retirar uma amostra de lodo do fundo da coluna para determinar o teor de sólidos totais de cada 
amostra de acordo com o método descrito em APHA, AWWA, WEF (2005). 


-Retirar uma amostra do sobrenadante e determinar cor aparente, turbidez e pH. 


-Esvaziar a coluna de adensamento e lavar com água o conjunto agitador e coluna. 


Nesta etapa você vai precisar: 

1 coluna de adensamento (já descrita neste item), 1 cápsula de porcelana de 100 mL, 1 dessecador 
completo, 1 estufa, 1 balde de plástico graduado com capacidade pra 5 L, 1 proveta de vidro de 100 mL, 1 
pisceta com água destilada, 3 béquers de 50 mL, 1 turbidímetro, 1 colorímetro, 1 medidor de pH e 1 
balança de precisão com 4 casas decimais. 


Não se esqueça de estar adequadamente vestido com calças compridas, sapatos fechados, jaleco, luvas e 
máscara. Também fiquem atentos para fazer registros fotográficos em todas as etapas. 
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9.4. TRATAMENTO DOS DADOS OBTIDOS LABORATÓRIO 


Todos os dados devem ser organizados em planilha eletrônica Excel ou equivalente e trabalhados para a 
determinação da velocidade de sedimentação e a melhor condição testada. 


= Preencher a Tabela 2 com os valores de TST obtidos no fundo da coluna de adensamento em cada 
condição testada. 


Tabela 2- Teor de Sólidos Totais (TST) no fundo da coluna de adensamento para cada condição testada. 


Sem polímero Com polímero 
T.S.T (%) E a T.S.T (%) Dosagem de polímero (g/kg) T.S.T (%) 
1 4 
n Preencher a Tabela 3 com os valores de cor aparente, Turbidez e pH obtidos em cada condição 
testada. 
Tabela 3 - Cor aparente, turbidez e pH para cada condição testada. 
Ensaio Cor aparente (UC) Turbidez (UT) pH 
Sem polímero 

Com polímero 1g/kg 

Com polímero 2g/kg 
n Gerar gráfico de altura de interface com linha de tendência ao longo do tempo para cada condição 
testada. 
n Calcular, para cada curva gerada, a velocidade de sedimentação em zona (vsz) pela tangente do 


trecho inicial retilíneo. Preencher a Tabela 4. 


Tabela 4 - Velocidade de sedimentação em zona (vsz) para cada condição testada 


Sem polímero | Com polímero 
Dosagem de polímero Dosagem de polímero 
Vsz (cm/s Vsz (cm/s Vsz (cm/s 
o (8/48) Sa (8/48) o 
1 4 
= Calcular os coeficientes A e B e estabelecer a equação da VSZ para cada condição testada. 
u Comparar os resultados obtidos visando explicar e definir a melhor condição testada. 
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10.1. OBJETIVOS INSTRUCIONAIS 


Ao final da aula prática de análise de oxigênio dissolvido em amostras de água o aluno estará apto a: 


1. Compreender a importância do oxigênio dissolvido como indicador de qualidade de ambientes 
aquáticos; 

2. Entender as etapas no desenvolvimento da análise; 

3. Realizar a análise de forma adequada; 

4. Calcular os resultados da quantidade de oxigênio dissolvido presente na amostra; 

5. Interpretar os resultados obtidos, comparando com as legislações ambientais vigentes, dentro das 


especificidades e objetivos de cada trabalho. 


10.2. A IMPORTÂNCIA DO OXIGÊNIO DISSOLVIDO (OD) COMO INDICADOR DE QUALIDADE 


O oxigênio dissolvido, dentre outros parâmetros ambientais, é fundamental na manutenção do equilíbrio 
dos ecossistemas aquáticos e na qualidade da água de mananciais. Muito utilizado em estudos ambientais 
é indicador de impactos relacionados aos lançamentos de efluentes, sendo um dos principais parâmetros 
de caracterização dos efeitos da poluição. 


É fundamental nos estudos de autodepuração natural de corpos d'água e no monitoramento e controle 
operacional em estações de tratamento de esgoto. Pelo fato de fornecer uma boa indicação da condição de 
ambientes aquáticos é utilizado com frequência no planejamento e gestão de recursos hídricos. 


O oxigênio pode ser inserido de ambientes aquáticos de várias formas, processos naturais e artificiais 
podem atuar para essa inserção de oxigênio no ambiente aquático. As plantas aquáticas e algas, durante os 
processos de fotossíntese, liberam e contribuem com oxigênio para os corpos d'água. A Atmosfera também 
contribui para aeração em corpos d'água, por difusão. Em mananciais superficiais com quedas d'águas 
e/ou corredeiras e intensa movimentação da água, a aeração é potencializada. Existem ainda, processos 
artificiais que são utilizados, principalmente, em processos de tratamento de efluentes. A depleção de 
oxigênio pode ocorrer a partir da oxidação de íons metálicos, da degradação de matéria orgânica, e pela 
elevação da temperatura de corpos hídricos. Esse rebaixamento da concentração do oxigênio em 
ambientes aquáticos é prejudicial à manutenção dos ecossistemas aquáticos. 


A resolução 357/05 da legislação do Conselho Nacional do Meio Ambiente - CONAMA estabelece o valor 
mínimo de 5mg/L (águas classe 2, quando não se te sem o enquadramento do corpo hídrico) como valor 
seguro para equilíbrio e manutenção de ambientes aquáticos. 


10.3. ETAPAS DO DESENVOLVIMENTO DA ANÁLISE DE OXIGÊNIO DISSOLVIDO (OD); 


A análise de oxigênio dissolvido pelo método químico, também conhecido como método de Winkler 
modificado ou iodométrico, é realizada em algumas etapas até sua finalização completa. As etapas 
seguidas no método são: 


Coleta da amostra 
Fixação da amostra; 
Liberação da amostra; 
Titulação da amostra; 


é vin 


A coleta de amostras para análise de oxigênio dissolvido pode ser realizada utilizando-se um batiscafo 
ou uma garrafa hidrológica. O batiscafo (Figura 1) é indicado para realizar coletas superficiais ou 
subsuperficiais até 30 cm da coluna d'água. 
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Figura 1 - Batiscafo para coleta de amostras superficiais 


Fonte: CETESB 2011. 


A água entra através de um tubo localizado na parte superior da tampa do equipamento e atinge o interior 
do frasco, fazendo com que o ar contido saia por um orifício lateral à medida que ele vai sendo preenchido 
com água (CETESB, 2011). Para coleta de amostras superficiais e em diferentes profundidades utiliza-se a 
garrafa hidrológica. Existem modelos diferentes e capacidades variadas de garrafas hidrológicas. Os tipos 
mais utilizados são van Dorn e Niskin (Figura 2). 


Figura 2 - Garrafas hidrológicas de coleta de amostra de água em superfície e em profundidade; A- garrafa 
de Niskin; B- garrafa de van Dorn 


cabo 
mensageiro 


Hberador 


Fonte: Adaptado de catálogo Hydro-Bios 


Para realizar a coleta de amostras a garrafa é mergulhada com ambas as extremidades abertas, ao atingir a 
profundidade desejada, o mensageiro (Figura 3) é solto e fecha o amostrador hermeticamente. 
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Figura 3 - Modelos de mensageiros de garrafas hidrológicas 


Fotos: Carlos Jesus Brandão /CETESB). 


Na amostragem utilizando a garrafa hidrológica, a amostra coletada deve ser distribuída para o frasco de 
DBO (frasco de vidro borossilicato com tampa esmerilhada e selo d'água) de capacidade de 300 ml (Figura 
4), de forma cuidadosa, a fim de evitar a aeração da amostra no momento da transferência. É importante 
que a amostra coletada não contenha bolhas de ar no seu interior, o que pode causar inferência e levar a 
um falso resultado. 


Figura 4 - Frasco de DBO com capacidade de 300 ml 


Nesse tipo de procedimento, é necessário introduzir a mangueira de saída de água da garrafa hidrológica 
até o fundo do frasco DBO e abrir o bico de drenagem devagar, enchendo o frasco cuidadosamente até 
extravasar a uma quantidade de água e só então tampar o frasco observando a ausência de bolhas de ar no 
interior do frasco. A incidência de luz também pode alterar as características iniciais da amostra e por isso 
é importante que após a coleta da amostra, o frasco seja mantido protegido da luz. 


A fixação do oxigênio dissolvido na amostra é realizada com a adição de 2 (dois) reagentes, o sulfato 
manganoso (MnSO4)(R1) e a solução álcali-iodeto-azida (R>), que é uma solução de hidróxido de sódio 
(NaOH), iodeto de sódio (Nal), e azida sódica (NaN3). Os reagentes Rj e R; são adicionados imediatamente 
após a coleta, a fim de evitar que ocorram alterações na concentração do oxigênio presente na amostra. 
Nessa etapa ocorre a formação de óxido manganês (Equação 1). Com intensa floculação da amostra 
(Figura 5). 


Mn (0H2) + 402 > Mn0O2 + H20 (1) 
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Figura 5 - Floculação na amostra causada pela adição dos reagentes R; e Ro. 


À azida sódica foi incorporada a análise para remover os nitritos presentes que são interferentes na 
reação, conforme equação 2,3; 


NaNs + H+ > HNs + Na* (2) 
HNs + NO2- + Ht 5 No + N20 + H50 (3) 


A liberação da amostra para a titulação (última etapa da análise) é realizada pela adição de ácido sulfúrico 
concentrado (Equação 4), que quebra os flocos formados e causa o aparecimento de uma coloração 
amarelada, que é proporcional a concentração de oxigênio dissolvido contido na amostra. 


MnO; + 21 + 4H+ > Mn*? + [; + 2H20 (4) 


Nesta etapa da análise o íon iodeto anteriormente formado, é oxidado a iodo molecular em quantidade 
proporcional ao óxido de manganês presente na amostra, que é também proporcional à concentração de 
oxigênio dissolvido presente. 


A última etapa da análise é a titulação do iodo liberado na amostra, com o reagente tiossulfato de sódio 
(Na2S406) (Equação 5). Esta etapa é também conhecida como iodometria. 


2Naz8203 + I; > Naz8406 + 2Nal + 10H20 (5) 
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A amostra é titulada com Tiossulfato de sódio (NazS203) 0,0125 N. Uma solução de amido é utilizada como 
um indicador de viragem nessa titulação. O ponto final da titulação é a viragem da coloração azul para 
incolor. 


10.4. PRÁTICA EM LABORATÓRIO 


Nossa aula prática acontecerá no laboratório ou em campo (com adequação de local para as titulações), 
com medições reais e na sequência de etapas do método e ao final da titulação o cálculo do resultado com 
possibilidade de comparação com os valores permitidos pela legislação ambiental vigente para o 
parâmetro analisado. 


É interessante ressaltar que a titulação para a determinação do oxigênio dissolvido em amostras de água 
tem um prazo definido para ser realizada, que é de até 8h após a etapa de fixação (etapa que preserva a 
amostra de alterações de suas características iniciais) (CETESB, 2011). E é por esse motivo, que em 
algumas situações é necessário realizar a titulação em campo. Nesse sentido, é possível levar até o local de 
coleta o material necessário a realização da análise, organizando uma bancada com as vidrarias e 
reagentes necessários para a análise com os mesmos cuidados, atenção, organização como no laboratório. 


10.4.1. SELEÇÃO DO LOCAL 


O local escolhido deve atender aos critérios apresentados em sala de aula. Para determinação da 
quantidade de oxigênio presente em um corpo d'água é necessário observar in loco um local onde é 
possível realizar a coleta de forma segura e com representatividade; observar dentro do objetivo do 
estudo, se as amostras são de águas superficiais ou subterrâneas e utilizar o equipamento indicado para a 
realização da coleta. Para executar as etapas já mencionadas anteriormente de forma adequada é 
importante organizar a bancada com todo material necessário ao desenvolvimento da metodologia. Na 
etapa de laboratório, a utilização de EPI (Equipamentos de Proteção Individual) é indicada, além de 
vestimenta adequada, como jalecos de mangas compridas, calça jeans, sapato fechado. Após o 
procedimento de coleta, acondicionamento e preservação das amostras segundo Guia de coleta e 
preservação de amostras de água (CETESB, 2011) separar o material necessário para análise. 


10.4.2 MATERIAIS/ REAGENTES; VIDRARIAS UTILIZADOS NA ANÁLISE 
10.4.2.1 REAGENTE UTILIZADOS NA ANÁLISE 

Solução de sulfato de Manganês; 

Solução alcalina de iodeto e azida sódica; 

Ácido sulfúrico concentrado; 

Solução indicadora de amido 1%; 


Solução padronizada de tiossulfato de sódio (0,0125N); 


10.4.2.2 MATERIAIS E EQUIPAMENTOS 
Batiscafo; 

Garrafa hidrológica; 

Frascos de DBO de capacidade de 300 ml; 
Buretas de 25 ml; 

Proveta de 100 ml; 

Pipeta graduada de 10 ml; 

Erlenmeyer de 250 ml; 
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10.4.2.3 PROCEDIMENTO ANALÍTICO 
1. Encher o frasco de OD com a amostra coletada, deixar extravasar um pouco (evitar aeração); 
2. Tampar o frasco tendo o cuidado de não deixar bolhas de ar no interior da amostra; 
3. Adicionar nesta ordem: 
Iml de MnSO,, 
Iml de Iodeto alcalino; 
Tampar e agitar por inversão por 3 vezes; 
Deixar decantar por cerca de 3 minutos; 
Iml de H;SO, concentrado. 


4. Tampar e agitar por inversão (a coloração amarelada indica a presença de oxigênio na amostra (Figura 
6), a presença de um precipitado branco, é indicativo da ausência de oxigênio dissolvido na amostra 
(Figura 7); 


Figura 6 - Presença da coloração amarelada indica a presença de oxigênio na amostra de água 


Fonte: Autora 


Figura 7 - Precipitado branco indicando ausência de oxigênio na amostra de água 


Fonte: Autora 


5. Passar 100 ml da solução que indica a presença de OD para um erlenmeyer de 250 ml (esta etapa é 
realizada em triplicata); 


6. Titular com NazS;03 (0,0125N) até que a solução fique amarelo pálido; 
7. Adicionar então 1 ml de indicador amido à solução; 


8. Continuar a titulação com o NazS,03 até ponto final de viragem de azul para incolor; anotar o volume 
total de NazS203 gasto na titulação. 
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10.5 CÁLCULO DOS RESULTADOS APÓS A FINALIZAÇÃO DA ANÁLISE; 


A expressão dos resultados para a análise de oxigênio dissolvido é dada em mg02/L e é obtido pela 
fórmula abaixo; 


Vix N(Na,8,05)x fx 8000 


O,(mg/L) = A 


Onde: 
V; = volume de Tiossulfato (NazS203) gasto na titulação; 
N = concentração da solução de Tiossulfato (Naz8203-0,0125N); 
F = fator de correção (obtido na padronização do Tiossulfato); 


V> = volume da amostra titulada. 


10.6 INTERPRETAÇÃO DOS RESULTADOS OBTIDOS, APLICANDO AS LEGISLAÇÕES 
AMBIENTAIS VIGENTES 


As resoluções que estabelecem valores máximos permitidos (VMP) para parâmetros ambientais em 
mananciais superficiais, são as resoluções do Conselho Nacional de Meio Ambiente - CONAMA. Em se 
tratando de água doces, classe 2, aplica-se o valor de 5mg/L de Oz, como valor mínimo necessário a 
manutenção dos ecossistemas aquáticos (CONAMA 357/05). Assim, ambientes aquáticos com valores 
inferiores a 5mg/L de O, podem indicar aporte de matéria orgânica no corpo hídrico, o que provoca o 
consumo de oxigênio, causando desequilíbrio no ecossistema. É importante ressaltar que uma única 
análise não é suficiente para diagnosticar impactos em determinado corpo hídrico. O monitoramento 
ambiental, com uma malha amostral bem definida e em períodos que seja possível verificar a variação 
sazonal do parâmetro ambiental estudado é importante para se ter uma avaliação mais segura e com isso 
interpretações mais corretas do ambiente estudado. 
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